Quantum Spin Nematic States in Bose-Einstein Condensates 
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We review some recent results on discrete symmetries and topological order in spinor Bose- 
Einstein condensates (BECs) of ^^Na. For spin one bosons with two-body scatterings dominated 
by a total spin equal to two channel, the BECs are in quantum spin nematic states at a low density 
limit. We study spin correlations in condensates at different limits and analyze hidden symmetries 
using a non-perturbative approach developed recently. We further more investigate the influence of 
hidden Z2 symmetries and U{1) quantum orders in symmetry partially retored states, particularly 
the effects on topological excitations. 



I. INTRODUCTION 



BECs of alkali atoms first created and studied experi- 
mentally in 1995-1997 revived theorists' interest on this 
subject [0J^,^. Because the Zeemann coupling in gen- 
eral provides stronger confining potentials for atoms com- 
pared with other trapping methods, and a rf evapora- 
tive cooling is needed to achieve BECs, the first gener- 
ation of BECs was realized in magnetic traps in labora- 
tories. In those early experiments, spins of atoms were 
completely polarized and BECs with single component 
were achieved, similar to the BEC of '^He studied a few 
decades ago. Apart from the complication of the con- 
fining potentials, the theoretical description of BECs of 
alkali atoms with positive scattering lengths in magnetic 
traps can be carried out in the Gross-Pitaveskii approach 
The activities on this subject are enormous and I 
should refer to an excellent review by Dalfovo et al. ||]. 

The physical properties of single component BECs of 
alkali atoms depend on two length scales. The inter- 
atomic distance, which is about a few hundred nano- 
meters in most of experiments, characterizes the kinetic 
energy in the problem. The BEC transition temperature 
which is determined by this length scale is about a few 
micro Kelvins or lower in experiments. The upper limit of 
the atomic density in BECs and the transition tempera- 
ture are imposed by three-body recombination processes 
which lead to trap losses in experiments 0. The sec- 
ond length scale is the two-body scattering length, which 
is about a few nanometers for alkali atoms. It charac- 
terizes the interactions between atoms and therefore the 
chemical potentials. Given the density of atoms and the 
scattering lengths in most of experiments, the chemical 
potentials are typically of the order of a few hundreds 
nano Kelvins. Amazingly, the scattering length depends 
on magnetic fields and the variation in scattering lengths 
is particularly substantial when the atoms are close to 
the Feshbach resonance 1^,^. This remarkable property 
of alkali atoms was observed in Na by Inouye et al. , 
in Rb by Courteille et al. ||ll|, Roberts et al. ||l2| and 
recently Cornish et al. p4[ , in Cs by Vuletic et al. [ p^ . 
Collapsing of BECs of "Li atoms with negative scatter- 



ing lengths was studied by Sackett et al. Dynamics 
of collapsing was investigated theoretically in and 
remains to be fully understood. 

Correspondingly, there are at least two classes phe- 
nomena investigated very extensively by experimental- 
ists. The first one is connected with quantum coherence 
of BECs. Quantum interference between two condensates 
(both spatial and temporal) and the creation of vortices 
belong to this category |l3,|l|,|l|,|o|,|ll . For instance, 
vortices were first created using a combination of inter- 
nal transition and external rotation by Matthews [|l9| , as 
theoretically suggested by Williams and Holland in 1999. 
Arrays of vortices were found by Madison et al. who 
stirred BECs by rotating the trap. Superfluidity of BECs 
has been examined by Raman et al. who demon- 
strated frictionless flow and by Marago et al. who 
observed beautiful undamped scissors mode suggested by 
Guery-Odelin and Stringari in 1999 p^ . 

The second class is related with scatterings between 
atoms, or the collective behavior of BECs. The AC 
Josephson effect of two internal states with different 
scattering lengths and zero sound excitations at coher- 
ent limit fall into this category j2^j2^,^,^. The sec- 
ond class phenomena are of particular interest from con- 
densed matter physics points of view and let me describe 
two interesting experiments in some details. In the ex- 
periment in P5[ , Hall et al. first prepared ^'^Rb in a 
hyperfine spin state \F = l,mp = —1 > in magnetic 
traps. By applying external microwaves or so called ^- 
pulses, they successfully transformed about fifty percent 
of atoms from |1,— 1 > state into |2,— 2 > state. The 
chemical potentials of these two condensates are differ- 
ent because of different scattering lengths. After certain 
waiting time T, a fraction of the condensate at |2, —2 > 
state is transformed back to |1, — 1 > state by another ^- 
pulse. The interference between the original condensate 
and the newly converted condensate depends on the ac- 
cumulated phases of two condensates after time T, which 
are different because of two different chemical potentials 
fii = Anaip/M for condensates with a density p, a scat- 
tering length ai and an atomic mass Af . The temporal 
dependence of the interference observed by Hall et al. 
does exhibit periodical oscillations at ms scale, expected 
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by a simple estimate of the chemical potentials. 

The observation of zero sound excitations in '^^Na 
BECs was carried out by Stamper-Kurn et al. in 
1998. After a temporal magnetic field was applied to dis- 
turb the BEC, a breathing of condensate was observed 
at a temperature much lower than the BEC transition 
temperature. Snap shots taken at 5ms apart indicate os- 
cillations in the size of the BEC with a period of a few ms, 
corresponding to the time interval for a sound like excita- 
tion to travel from one end of the BEC to the other one. 
The sound velocity, being proportional to the square root 
of the compressibility, depends on the scattering lengths 
and appears to be consistent with the frequency of the 
oscillations. These oscillations were also shown quickly 
damped at the temperature close to the transition point, 
caused by scatterings with a thermal cloud. 

BECs with multiple components were created and 
studied a few years later [p9pc|j3l|] . Particularly, 
Stamper-Kurn et al. successfully used optical beams to 
trap '^^Na, which confines alkali atoms independent of 
spins and liberates spin degrees of freedom Q . For alkali 
atoms '^^Na or with nuclear spin / = 3/2 and elec- 

trons at s orbits, the energy splitting between hyperfine 
multiplets is of order lOOmfc. At temperatures as low as 
a few hundred nano Kelvin, '^^Na and atoms can be 
considered as simple bosons with a hyperfine spin F = I. 
Liberation of spins in BECs in an optical environment 
indeed leads to a new dimension of studying BECs, espe- 
cially spin correlated condensates. This further inspired 
many theoretical activities on spinor Bose-Einstein con- 
densates mjUg. The spin correlated states in ^^TVa 
were afterwards studied experimentally by Stenger et al. 

The new complication in BECs in optical traps surely 
arises from the extra internal degree of hyperfine spins. 
When Nmp{ mp — 0, ±1) atoms condense at each of the 
hyperfine states \F^mF > but all condense at the low- 
est orbital state, the ground state of N{= "YlmF -^rnp) 
noninteracting spin F ^ \ bosons has a I?s-fold spin 
degeneracy. 
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by contrast to its magnetically trapped cousins. Pre- 
sumably, hyperfine spin-dependent scatterings lift the de- 
generacy and lead to spin correlated states. Optically 
trapped BECs therefore set up a platform for the study 
of the quantum magnetism in many boson systems and 
add a new twist to already extremely fascinating systems. 

For ultracold atoms in BECs, the spin-dependent two 
-body interaction is characterized by ?72(ri, y^^) — 5{yi — 
r2)[co+C2FiF2], as suggested in an early paper |3^. Here 
Co = (50 + 2g2)/3, C2 = (32 - 5o)/3; qf = ATTh^ap/M, M 
is the mass of the atom and a^? is the s-wave scattering 
length in the total spin F channel. Thus, spin corre- 
lation in BECs is determined by C2. For ^"^Rb studied 



experimentally by the JILA group ||2^ , 52 < ffo or C2 < 
and the scattering is dominated by the total spin F = 
channel. In the ground state, all spins of atoms prefer 
to align in a certain direction and have a maxima mag- 
netization [32[. For ^'^Na studied experimentally by the 
MIT group ]30[|, the scattering between '^■^Na atoms is 
dominated by the total spin F — 2 channel, i.e. (72 > go- 
The scatterings between ^"^iVa atoms thus lead to "an- 
tiferromagnetic" spin correlation. In this article, I will 
fully explore the properties of the BEC of '^^Na. 

Efforts have been made to understand the ground state 
properties, exact excitation spectra, collective modes and 
topological excitations. Theoretically, issues of spin or- 
dering in the spinor BECs are of particular interest and 
were treated by Ho, and also Ohmi and Machinda in 
1998. In both works, the spin ordering was discussed in 
a mean field approximation, in the context of the Gross- 
Pitaveskii equation. Ho also obtained two branches spin 
wave excitations by a Gaussian expansion. 

The first attempt to understand spin correlation be- 
yond the mean field limit was made by Law, Pu and 
Bigelow [|4|. Law et al. observed that the micro- 
scopic Hamiltonian of N spin-1 Bosons interacting via 
a spin-dependent two-body potential given above com- 
mutes with the total spin of the system; the spin is a 
good quantum number and the exact eigenstates should 
be eigenstates of the total spin operator. When C2 > 0, 
the ground state should be a spin singlet and rotation- 
ally invariant. This was explicitly demonstrated in a zero 
mode approximation, where the finite momentum sector 
of the Hilbert space was neglected Q . Using a four- wave 
mixing method. Law et al. obtained the exact excitation 
spectrum in a OD limit, a limit where all spatial fluctu- 
ations are neglected. These results in OD case were also 
derived by Castin and Herzog |^5|, using a geometrical 
description very close to the one I am going to introduce 
in this paper. Law et al.' s results, though are barely 
relevant to the particular experiment done by Stenger et 
al., given the number of atoms, time scales during which 
the experiment was performed and small biased magnetic 
fields in the experiment, nevertheless inexplicitly reveal 
the role of symmetry restoring processes and the nature 
of quantum fluctuations of the spin order in the problem. 
The even-odd effect in the excitation spectrum observed 
by Law et al. further implies a discrete symmetry in the 
problem which I believe has a profound origin. 

Ho and Yip later found that the singlet states obtained 
by Law et al. are rather sensitive to an external magnetic 
field 1^^. In a thermal dynamical limit, an infinites- 
imal magnetic field will stablize the symmetry broken 
state and the rotationally invariant state is irrelevant as 
far as experiments with a small magnetic field are con- 
cerned. This feature of '^^Na was later realized to re- 
late with the level spacing of the low lying excitations 
P^,p8|. For the OD dimensional BEC, the energy of low 
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lying excitations scales down as the number of atoms and 
the quasi-degeneracy makes the system ultrasensitive to 
magnetic fields. A symmetry broken state, which can be 
constructed out of low lying excitations has an energy in- 
versely proportional to the number of atoms involved and 
can be pinned by an infinitesimal field when the number 
of atoms is infinity. 

Topological excitations in the spinor BECs of ^"^iVa 
have been paid a special tribute to by a few groups. Lin- 
ear defects were studied by Leonhart and Volovik [ p9[ , 
who pointed out the existence of Alice strings in the 
BECs of "^^Na. Textures of the Skyrmion type were 
pointed out in [^2| ; and most recently, point defects were 
examined [0. Except the work of Leonhart and Volovik, 
most of the works on this subject were based on Ho's 
original identification of the internal space, which missed 
the important discrete symmetry of Z2 type. Therefore, 
many new aspects of topological defects studied in this 
paper have never been investigated. For instance, in a re- 
cent work pol , the authors applied the nonlinear a model 
derived in |37[|3q ] to study static and dynamical proper- 
ties of hedgehogs. Unfortunately, the authors of that 
work failed to observe the correct internal space. Many 
aspects of the point defects in the BECs of ^^Na, particu- 
larly those under the influence of the entanglement of the 
phase and spin orders including the homotopical indis- 
tinguishability between ± hedgehogs, between hedgehogs 
and closed loop 7r-disclinations were previously discussed 
prf . I will review those results in this article. 

The pursue of having a nonperturbative approach be- 
yond a zero dimension approximation was first made in 
That approach allows one to identify that the 
low energy spin dynamics in the BECs of spin one atoms 
with antiferromagnetic interactions belongs to the same 
universality class of the nonlinear-cr model(NL(TM), a rel- 
atively well known subject, thanks to many studies in the 
context of field theory. The NLcrM also became known 
to many condensed matter physicists working with anti- 
ferromagnetic systems after a paper in 1983 ||4l|] . 

In ||37| , |38|| , the spin dynamics in the problem of N spin 
one bosons interacting via a spin-dependent two-body 
potential with C2 > was mapped into an o{n) nonlin- 
ear sigma model (NLcrM). n = 3 at the zero magnetic 
field limit and n = 2 in the presence of a weak magnetic 
field. With the help of this mapping, the many-body 
aspect of '^^Na becomes rather explicit and many prop- 
erties including spin correlations in Id can be studied. 
Furthermore, energetics and dynamics of topological ex- 
citations in BECs can therefore be easily analyzed in this 
approach. 

At a low density limit, it is also found that the ground 
state has a similar discrete symmetry as that of a uniaxial 
nematic phase of liquid crystals, except a 7r-phase factor 
as discussed extensively in p^ , |38[ | and for this reason the 
condensate is at a quantum spin nematic stSLte{QSNS). 
The ground state degeneracy space is [U{1) x S^]/Z2, 



with n depending on a magnetic field. The antiferro- 
magnetic scatterings between cold atoms lead to a Z2 
symmetry breaking. 

This discrete symmetry was later appreciated by Dem- 
ler et al. in the context of low lying excitations | ^ , ^ . 
Demler et al. noticed that the Z2 symmetry in a sym- 
metry broken state actually also appears in all low lying 
states. By imposing a proper projection in a functional 
integral approach, they showed that the full low energy 
physics in the phase and spin sector is given by Z2 gauge 
fields coupled with two matter fields characterized by an 
XY model and a NLctM. This work as well as the pre- 
vious work [ p7| , ^ open a new door for the studies of 
many-boson states. Many issues, such as fractionaliza- 
tion of topological excitations, fractionalization of quasi- 
particles and topological order were studied and investi- 
gated with the help of this mapping Qj45|. A number 
of new predictions have been made for the spinor BEC, 
most of which have never been studied experimentally. 
The mapping also provides a starting point for the future 
study of some strongly correlated fcrmionic systems. 

The purpose of this article is to provide a guideline and 
introduction to some ideas on the spinor BECs developed 
very recently in Through out the article, I will 

emphasis on the concepts of topology, at levels of both 
the mean field approximation, which is suitable for a dis- 
cussion on symmetry broken states and beyond, which is 
essential for the study of symmetry restored states. This 
point of view might be dramatically different from the 
conventional ones shared by atomic physicists. On the 
other hand, it seems to me that the nature has chosen to 
behave in such a unconventional, complicated way that 
I am left with no alternatives, or a simpler presentation. 
It also appears to my collaborators and me during this 
investigation that this is a subject full of surprises and 
new issues keep on poping up. At the time of writing, 
I present results which I believe to be relevant to BECs 
investigated in experiments. There are many other is- 
sues still remaining to be explored. Furthermore, some 
descriptions given here are qualitative and conceptional, 
so as to be communicated in a most efficient way. Like 
many other new subjects, a more quantitative aspect of 
some properties can be achieved only after basic notions 
are established and will be present in a series of separate 
works. At last, I should point out that many strongly 
correlated electronic systems are believed to share very 
similar topological properties and therefore the 

BECs might be systems where some of basic ideas es- 
tablished in strongly correlated electron systems can be 
tested experimentally. Understanding of the spinor BEC 
might as well shed some light on other Fermionic super- 
fluids such as p-wave triplet superconductors. 

So, I will review the nonperturbative approach re- 
cently formulated. Results obtained previously by dif- 
ferent groups with other methods can be reproduced in 
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this new representation. I will of course go beyond previ- 
ous framework and outline many novel results obtained 
in the last two years |^,||,|U4|g|,||l . In section II, a 
mapping from the original microscopic model of N spin 
one Bosons with C2 > into a NLaM at long wave length 
limit will be introduced. Furthermore, the corresponding 
N-body wave function is shown to be invariant under a 
180° rotation of a vector n and a U{1) gauge transforma- 
tion. At last, Z2 gauge fields are introduced in the model 
to enforce the discrete symmetry. 

In section III, symmetry broken states are examined. 
The internal order parameter space for the BEC is identi- 
fied as X 5^]/Z2(zero field limit). Linear, point defects 
and " particle like" textures are discussed with an empha- 
sis on the influence of the group and the discrete Z2 
symmetry. Particularly, a family of strings representing 
TT spin disclinations superimposed by superfluid vortices 
of half integer circulation are studied. The distinguisha- 
bilities of ± hedgehogs, both physical and homotopical, 
are analysised. The connection between closed loop tt- 
disclinations and hedgehogs is established. 

The effect of an external magnetic field is also re- 
viewed. The NLcrM in this case has two components 
instead of three at zero external field. The ground state 
is degenerate in a manifold [C/(l) x S'^]/Z2. 

In section IV, symmetry unbroken states are studied 
in the scheme of the NLcrM. The possibility of having 
ncmatic order-disorder phase transitions and the nature 
of spin correlation are addressed. The ground state of 
N spin one Bosons interacting via C2 > spin-dependent 
two-body potentials in Id limit is also present. In section 
V, the influence of the Z2 symmetry on the symmetry un- 
broken states is discussed. One of the most interesting 
consequencies of having Z2 symmetries is the possible 
fractionalization of topological excitations. Under the 
influence of the Z2 symmetry, it is shown that vortices 
can be fractionalizcd into half-vortices. It is also pointed 
out that instead of a condensate of individual atoms, one 
can have a singlet pair condensate, which supports free 
Z2 vortices. 

In section VI, the property of connection fields is anal- 
ysised and hidden topological order is identified in spin 
disordered BECs. The existence "topological order" is 
argued due to an "order from disorder" mechanism. 

In section VII, I discuss the symmetry restoring of a fi- 
nite size nematically ordered state. Quantum mechanical 
nature of the order parameters in finite size ^^Na BEC 
was addressed in connections of experimental probes, 
taking into account the metastability of the atomic gas. 
In section VIII, I make some remarks on the relation 
between the BECs considered here and other strongly 
correlated magnetic systems. 



II. GEOMETRIC DESCRIPTION 

The inconsistency between the mean field point of view 
taken in te2|j3^] and the spin singlet ground state ob- 
tained in ||34| , lies in the heart of the quantum symmetry 
restoring. If one starts with the symmetry broken solu- 
tion in the CPE approach, for instance, one can ask how 
the rotation symmetry is restored after all and when that 
happens. A more practical question would be, given the 
life time of the metastable BEC gaseous cloud studied in 
the experiment, is the rotationally invariant state rele- 
vant? This question motivates us to develop a nonper- 
turbative approach. The answer to this question requires 
a complete characterization of nonlinear spin dynamics. 

The second issue concerning us is the validity of the 
zero mode approximation p4| . Though the exact solu- 
tions in a zero mode approximation provide information 
of low lying excitations, they fail to yield spatial-temporal 
spin fluctuations, which are essential for the most general 
description of the ground state of N interacting F — 1 
Bosons. The finite momentum sector of Hilbert space is 
clearly critical for such a discussion. 

It is also worth commenting that the CPE approach 
with three components |^,^, on the other hand, does 
permit a gradient expansion around the symmetry bro- 
ken solution in principal. However, such a procedure is 
feasible only when the derivation of the true ground state 
from the uniform solution is small. Practically, the GPE 
approach works in a small vicinity of the symmetry bro- 
ken solution in the order parameter space and is far from 
being a general description. 

The third and perhaps the most important one is the 
description of nonlinear spin dynamics in the spinor BEC. 
The failure of the generalized GPE approach far away 
from its mean field solution is largely due to the fact 
that the collective variables chosen there are not suitable 
for a description of the nonlinear spin dynamics. For this 
reason, in the framework of the GPE approach, the inter- 
actions between low lying excitations (spin waves) which 
eventually disrupt the symmetry broken state are barely 
tractable. 

For instance, when the zero point rotation energy of 
each atom caused by two-body scatterings is high enough, 
long-range spin order predicted in the GPE approach is 
destroyed even at zero temperature (see details in section 
IV). In that case, spin correlations can be nontrivial and 
yet not characterized by the simple solutions of the GPE. 
In addition, the spin correlated state in a very anisotropic 
trap, or Id, is clearly beyond the validity of the GPE ap- 
proach. 

Efforts along this line had been carried out recently. 
To have a general approach of studying spin correlated 
states of spin one bosons with C2 > was the purpose 
of those works. It is reviewed in the rest of this section. 
I will demonstrate a geometrical approach which allows 
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us to go far beyond the three-component Gross-Pitaveskii 
approach and explore a variety of possible spin correlated 
states. 



A. NLcrM characterization of spin dynamics 

Following the discussions above, to study the spin cor- 
related BEG, it is important that a set of collective co- 
ordinates are correctly defined. Idealistically, these vari- 
ables should be chosen to a) completely characterize the 
nonlinear spin dynamics in the whole parameter space 
and b) to form a simplest representation for two-body 
spin dependent scatterings. For this purpose, we use 
two vectors (f2i,r22) living on two unit spheres to de- 
scribe the nonlinear spin dynamics in the presence of 
spin-dependent scatterings. We are going to map the N 
interacting spin one Boson problem into an 0(3) NLctM 
self-consistently. 

We make two assumptions to proceed further. First, 
there exists a class of spatially varying N-body wave 
fimctions, out of which the N interacting spin-1 Bosons 
ground state and the low lying excitations can be con- 
structed. The whole class is characterized by a few col- 
lective coordinates( r2i.2 introduced below). 

Second, the energy of this class N-body wave functions 
under consideration can be expressed in terms of a finite 
number of collective variables. Quantization of the collec- 
tive modes yields elementary low lying excitations, which 
can interact with each other strongly. The exact N-body 
wave function of the ground state is characterized by the 
zero point motions of collective variables and their multi- 
couplings. Practically, in most limits, the ground state is 
conveniently characterized by the correlation functions. 
The low energy spin excitations are of a collective na- 
ture; the excitations of individual spins are neglected in 
the leading order approximation 0{1/N). 

The correctness of this approach is justified by the fol- 
lowing measure: The resultant theory which is written 
in terms of collective variables instead of N-body wave 
functions of individual atoms should produce the correct 
nonlinear spin dynamics, such as Euler's equation of mo- 
tion, conservation laws, at the long wave length limit. 
The theory is self-consistent if this measure is satisfied. 
The derivation of the mapping was obtained and sum- 
marized in | |37| , ^ . 

To describe the spin correlated BEG, it is most conve- 
nient to introduce Weyl representation of SU{2) involv- 
ing polynomials of a unit vector (m,w) |4^,^^. Each unit 
vector is represented by a point 17 on a sphere with polar 
coordinates (0,0); namely 



u = exp(i-)cos-,w = exp(-i-)sin-. (2) 
The corresponding hyperfine spin operators are 



F+^4,F-.4; 
ov ou 

^ 2 du dv' 



(3) 



The scalar product between two wave functions g and / 
is defined as / g*{u,v)f{u,v)(Kl/AT:. (We reserve for 
the spin rotation discussed below.) 

Under spin rotations = exp(zF2Xi/2) exp(iFy0i/2) 
exp(z_F'z(/)i/2), u and v transform into 

exp(i y) i Y exp(-« y)u + sm — exp(« y)?;) 
^'(^^2,X2) = 

exp(-« yH ~ si'^ y exp(-iy )m -I- cos — eyi-p[i — )v) 

(4) 

where r2i,2 — (6*1,2,01,2)- The following identities hold 



«*(f2i)F«(Jli)^ = if2i, 
47r 6 

j z;*(ni)F«(fii)^ - -if2i 



(5) 



Spin one wave functions are polynomials of degree 
two in u and v. \/3u^, ^/Quv, ^fiv^ correspond to 
m — 1,0,-1 states. All F — 1 states can also be ex- 
pressed in term of ^/Qu[Q,\)v{Q,'2) with r2i^2 properly 
chosen. The examples are given in Fig.l. The compari- 
son between the polynomial representation and the usual 
vector representation is given in the Appendix A. 

The Hamiltonian for spin one bosons can be written 

as 



a a, 13 

+ Y,^-^9tJ^BH. (6) 



The second term is the hyperfine spin-dependent inter- 
action and the last term is the coupling with an external 
magnetic field H — Hsz ; g is the g- factor of bosons and 
/is is the Bohr magneton. 

The wave function for N spin one Bosons can generally 
be written as 

*({ra}) = Sna^i„,N^N^{ra)VQUa{flla{ra))Va{^2a{ra))- 

(7) 

S is to symmetrize the wave function; labels one- 
particle states. The phase [xi — X2]/2 corresponds to 
a gauge transformation in the complex field $jv^ in- 
troduced above. Without losing generalities, we set 
Xi = X2 = 0. 
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By taking f2ia,2a = ^lA^) and $Ar„(r) = ^(r)/VN 
(<i>(r) is a complex scalar field), one chooses J7i^2,*&(r) 
as the collective variables of the N interacting spin one 
Bosons. One obtains a Hamiltonian for the spin and su- 
pcrfluid sectors 

H = Hs + He + Hex, 

'^s = lj dr[^{Wn{r))H*<f + 4c2L^(rW<i>\% 

'^c = lj dr[^|V$(r)|2+4co|<i>*(r)$(r)p], 

Hex = gtisH J dr$*(r)$(r)e, • L(r), (8) 

Here n(r) = {Qi + ^2)72, L(r) = (^i - fl2)/2. Sub- 
scripts s, c and ex label spin, charge and external fields. 
In the spin sector, the first term is the energy cost of 
twisting n in space, representing the spin stiffness; the 
second term is the effective " rotation" energy in the pres- 
ence of a finite spin moment L. 

A local spin-phase coupling term, Hsc, which repre- 
sents the superflow in the BECs due to the Berry's phases 



under spin rotation, does exist |37|. This term, however, 
is linear in terms of L and quadratic in terms of spatial 
gradient and is much smaller than Hs , He at the long 
wave length limit. Moreover, Hse vanishes in a state 
where L is zero. Based on these observations, it was 
pointed out that the local coupling is negligible p8|] . 

In the most interesting limit, we can introduce $(r) = 
■\/ p{r) exp(ix(r)); the local spin density is l(r) = 
L(r)p(r). n and 1 satisfy the constraint 



n(r) • l(r) = 0. 



(9) 



p and (p, n(r) and l(r) obey the following commutation 
relations, 

[p(r),x(r')] = iM(r-r'); 



[n,(r),n^(r')] = 0, 

[U(r),l0(r')] =z;ie"'''^V(r-r'). (10) 



a/37 



is an antisymmetric tensor. The second identity in 
Eq.|l^ is valid only when L per atom is much less than 
unity and n(r) can be considered as a classical "vector". 

The equation of motion for the superfluid and the spin 
sector can be derived as 

5a = 4cop-^VV+^(Vp)V, 
9tP = V.(^Vx), 
9,n = 4c2nx(l-^), 



4C2 



where the low frequency and long wavelength limits have 
been taken. J = [x, p~^]/ih. Eqjl^ determines the non- 
linear spin dynamics in the BECs in the long wave length 
limit. 

The first two equations indicate the finite compress- 
ibility or chemical potential in the presence of scatterings 
between atoms, and the usual current conservation law. 
The first equation suggests the AC Josephson effects in- 
volving two condensates in symmetry broken states. It 
also implies the noncommutative nature of the phase de- 
gree freedom x a-nd the Hamiltonian He for interacting 
atoms. This leads to the quantum fluctuations of the con- 
densate phase X and symmetry retoring in a finite con- 
densate. Linearizing these equations yields excitations of 
the nature of zeroth sounds, which involve a periodical 
compression of the density and a periodical modulation 
of the phase of the condensation. The sound velocity is 

Ve = y^4:CQp/M. 

The third equation shows that n precesses in the pres- 
ence of the spin density l(r) due to the two-body hyper- 
fine spin dependent scatterings. Again, it results from 
the noncommutative nature of the spin order n and Hs- 
The profound phenomenon of the quantum symmetry 
restoring of n in a finite size spinor BEC is driven by 
this property. In homogeneous BECs, the fourth equa- 
tion can be considered as the "conservation law" for the 
spin density l(r); the corresponding spin current density 
is jo, = pM^^ea/3rii^T]'^n0. p is the density of atoms. 
For the sake of simplicity, we will assume the system is 
homogeneous from now on. 

When C2 > 0, the "rotation" energy in Hs is positive. 
At the zero field limit, there exists a saddle point solution 
of Eq.H for the spin sector 



n(r) =no,l(r) = 0. 



(12) 



(11) 



no lies on the unit sphere. The solution is schematically 
shown in Fig.l. 

u(no)u(no) represents state (0,1,0)"^ with the quanti- 
zation axis pointing along ng. It is easy to confirm that 
this solution at 7? = corresponds to the "polar" state 
found in Ref. |^ . By expanding Eq.|ll| around the mean 
field solution, we obtain the spin waves with a sound 
like spectrum u; = ^ Ac2p/Mk, which can also be ob- 
tained in the GPE approach. The spin wave velocity is 
Vs = sjApcijM. 

However, only the vicinity of point no on the unit 
sphere is practically accessible in the representation em- 
ployed in Ref. Eq]|, on the other hand, is valid for 
any point Vl\ ~ on the sphere. Therefore, the ef- 
fective NLcrM derived below allows us to describe a spin 
correlated state far away from the one given in the GPE 
approach, as promised and provides at least a qualitative 
picture about the N-body ground state wave function 
well beyond the simple mean field solution. 
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(a) (b) (c) 



(d) (e) 




(f) 

Fig.l 1) States represented by two unit vectors ili and 
r22): a) u^; b) uv; c) v"^. 2) fli and ^2) in zero field d), 
and in an external field H e); 3) An example of micro- 
scopic wave functions of spin nematic states as indicated 
in f). In a)-f), arrows stand for J7i_2; shapes of shaded 
blobs are given by r = 3|u(J7)f in a polar coordinate 
(r, 9, (j)). Note in d) the orientation of n = (Jli +Q,2)/2 
is arbitrary while in e) n is perpendicular to the external 
magnetic field. 

The low frequency sector Lagrangian for the spin and 
the superfluid component can be derived by taking into 
account Eqs.| 11 

(13) 

We introduce r = as the imaginary time. The non- 
linearlity is imposed via the constraint |n'^| = 1 at a low 
frequency limit. Eq.^ is the main result of the map- 
ping and we keep terms which are of the lowest order in 
terms of dr and V. £c is taken in a Gaussian approxi- 
mation and should be replaced by a full Gross-Pitaveskii 
Lagrangian in general. We will be mostly interested in 



the physics in the spin sector and simplification in Cc 
doesn't affect conclusions in this paper. Cg in Eq.|l^ rep- 
resents an 0(3) NLcrM; a coupling between the spin and 
superfluid component due to the Berry's phase effect has 
been neglected p7| , ^ . 

The partition function of the Lagrangian can be con- 
sidered as the ground state wave function, with all the 
zero point fluctuations and their interactions included. 
The Lagrangian Cg has been been studied extensively 
before. Many fascinating properties of the spinor BEG 
therefore can be obtained by the mapping presented in 
this section. Depending on the ratio between two-body 
scattering lengths and interatomic distances, Eq.|l^ ad- 
mits solutions with very different spin correlations. The 
symmetry broken states will be discussed in the next sec- 
tion and symmetry unbroken states are to be addressed 
in section IV. 

B. Z2 symmetries and Z2 gauge fields 

In the previous subsection, we have shown that the 
low energy spin dynamics of the BEGs of spin one atoms 
with antiferromagnetic interactions is equivalent to that 
of the NLctM. These BEGs support two branches spin 
wave excitations. In this subsection, we are going to il- 
lustrate an extra discrete symmetry enforced by the an- 
tiferromagnetic interactions. This discrete symmetry ef- 
fectively couples the spin and phase degrees of freedom. 

Let us first look at the the ground state wave function 
shown in Eq.^ with ^1,2 — n. It is invariant under a 
global transformation n, x ^ — n, x -I- tt; i.e., 

«'(n, x) = *(-n, X + *(n) - {-lf-^[-n), (14) 

where x is the phase of the scalar field ^{x) intro- 
duced in Eq.0. In obtaining this symmetry, we notice 
u(n) = exp(i7r/2)u(— n), with 7r/2 from a phase of a spin- 
i particle under a 180" rotation. 

This is certainly a quite general feature of all low lying 
excitations in ^^Na with antiferromagnetic interactions. 
It can also be illustrated using the zero mode Hamilto- 
nian 

'Hz.^.=pc2 — , (15) 
L is the total spin of the BEG. And 
[na,n;3] = 0, 

[La,np] = ihe^'^'^n-^, [L„, L^] = ihe^'^'^L-^. (16) 

Eqs. p^ , p^ show that L is an angular moment operator 
defined on the unit sphere spanned by n. In a spherical 
coordinate where n = (sin 9 cos (f>, sin 9 sin <j>, cos 0), L is a 
differential operator with respect to 6*, </> |37j. The eigen- 
states of the zero dimension Hamiltonian are spherical 
harmonics of the collective variable n 
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*(n) = Yi: 



(17) 



The energy spectrum is given by — l{l + 1), which 
is identical to that in |Q. The symmetry of Bosonic 
wave functions imposes further constraints on parities of 
the states: Z = 0, 2, 4, ....N , if N is an even number; and 
I = 1,3,5, ... otherwise. Eq. |l^ shows the ground state 
is Z = or a spin singlet when N is an even number. For 
an odd N , the ground state has a threefold degeneracy 
with 1 = 1. The collective behavior of the ground state 
depends on the even-oddness of the number of atoms in 
the BEC. 





(a) 



(b) 



Fig. 2 The ground states for an even a) and odd b) num- 
ber of atoms in terms of the collective variable n. In a), 
the ground state is a non-degenerate s-orbit while in b), 
the ground state has a threefold degeneracy correspond- 
ing to one of px,Py,Pz orbits, following Eq.l7. 

Now we are ready to construct low lying wave pack- 
ets, taking into account the symmetry of Bosonic wave 
functions, 



*w(n) 



1=0,2. 
1=1,3. 



^/,mi1,m(n), N is even; 
Ai^rnYi^m{n), is odd. 



(18) 



One easily verifies that the wave function written in 
Eq.^ observes the same Z2 symmetry as a symmetry 
broken state. This is not surprising from the point of 
view that a symmetry broken state can be considered as 
a wave packet constructed out of the exact spectrum. 

The effect of the Z2 symmetry is manifestly dramatic 
on the dynamics of the spinor BEC. A full description of 
the BEC based on the Z2 gauge fields was given in pre- 
vious works, taking into account the entanglement of the 
order parameter space. So far, we haven't been able to 
impose this constraint of Eq.^ in the function integral 
in a continuous limit. However, by introducing an opti- 
cal lattice, we are able to enforce the Z2 symmetry in a 
mapping. A detail derivation was presented in a preprint 
by Demler et al. [^; a short account of results was given 
by Demler and Zhou |43j. I will illustrate the spirit of 
the derivation and quote the results here. 

Let me describe an optical lattice where the dynamical 
effects of the discrete symmetry can be conveniently dis- 
cussed [|0|j53j5|,|5|,||,|5p|]5i,|l|. Experimentally, fol- 
lowing Grynberg et al. pO|, a body-center cubic struc- 
ture can be created by having the following four laser 



beams interfere with each other: a circularly polarized 
light propagating along z direction and three linearly po- 
larized beams, with polarization plane coinciding with 
the xy plane propagating in directions perpendicular to 
the z axis. The lattice constant in an optical lattice is 
determined by the wavelength of the laser. In Gryberg 
et al.'s experiment, this is 514nm; the elasticity of the 
lattice is tunned by the laser intensity. So the model I 
am going to study is a realistic one from the experimen- 
tal point of view. On the other hand, all the phenomena 
we are going to discuss should occur in both single traps 
and optical lattices, though it is perhaps more practical 
to observe them in a lattice setting. The feasibility was 
discussed in Q. 

So we will have zero dimensional BECs with certain 
number of atoms living at each site described by Eqjl5|, 
subject to the constraint in Eq.[l8|. The atoms can also 
hop between sites in the presence of quantum tunneling. 
The Hamiltonian of an array of identical optical traps is 
given by 



i i'j 

-2 J n^Uj cos{xz - Xj) 



where 



d 

dXi 



(19) 



(20) 



and J w tN . The Hamiltonian given here also includes 
the phase degree of freedom to preserve the discrete sym- 
metry under a discrete gauge transformation. 

As indicated before, for each site characterized by the 
Hamiltonian in Eq.|l^, there are two sorts of elementary 
excitations: a)the excitations carrying no L but with iV; 
b)the excitations carrying no N but with L. All excita- 
tions are classified by two quantum numbers N and L as 
{n,l). The symmetry of the wave function observed by 
the low lying excitations is imposed by a constraint that 
the sum of n and I is an even integer. So only (n, 21 — n) 
are the physical excitations. For instance, a (0, 2) excita- 
tion corresponds to flip the spin of one of atoms without 
changing the number of the atoms; a (2,0) excitation is 
to add a singlet pair to the condensate. And at last, a 
(1, 1) excitation represents adding one atom of spin one 
to the BEC. 

To carry out a functional integral in the representa- 
tion of (ui, li) while preserving the symmetry of the wave 
function, we use a projection operator 



^ exp[i-(l-CT,)(T 

T,=±l 



(21) 



that introduces a new Ising variable Ui 
ical Hamiltonian therefore is 



±1. The phys- 
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(22) 



III. QUANTUM SPIN NEMATIC BECS 



Projecting an enlarged space (n, I) into {n, 2m — n) 
therefore results in a set of discrete variables which were 
identified as the Z2 gauge fields. The long wave length 
physics is characterized by the following action, 

5 = - ^ J^^,(Jrr'COSXrr' " ^ J^^, C0s{2Xrr') 
rr' rr' 

rr' rr' 

Here 

Qf = n>^ - - 5"^ (24) 
n 

is a nematic order parameter for the n-component unit 
vector n (n = 2 when n lies in the plane and n = 3 
when it can rotate in all three directions). The cou- 
pling constants are J^rif = Jtr±f — (<^5)~^j 
■K:t±{i,y....} - ^-n-ni and Jj^^^^^ ='j.';±f = 0, 
j'^'^'^t- - 1 = —eJ/A. e is the time unit introduced when 
we discretize the time; r = (i, r) and r + f = (i, t + e). 
And Xrr' = Xr — Xr' ■ Finally, 77 is a parameter in the 
Hubbard-Stratanovich transformation |^ , ^ . 

The site variables exp(ix) and n live in and 
respectively, representing the phase and spin degree of 
freedom at each site, arr' = ±1 are discrete linking vari- 
ables living in a Z2 space. In the 5^ sector, we recovered 
the NLcrM description of spin dynamics obtained in the 
previous sector. However, under the influence of the Z2 
symmetry, the spin-phase degrees of freedom are entan- 
gled, with interactions mediated by the Z2 gauge fields. 

It is obvious that the form of the action in Eq.|2^ is 
the simplest one consistent with the charge U{1), spin 
SO{3), and gauge Z2 symmetries of the model. Another 
term that is allowed by the Z2 symmetry and that is 
generated by integrating out the high energy degrees of 
freedom is the analogue of the Maxwell terms for the 
lattice gauge models 

Sa = -Kj2l['^rr' (25) 

□ □ 

where the summation goes over plaquettes in d + 1 di- 
mensional lattices. One can confirm that Eqs|2^j2^ are 
invariant under the following Z2 gauge transformation 

X ^ X + ^ ^, <^rr' ~> ~(Jrr'- (26) 

Note that the terms proportional to j2c,2s invari- 
ant under the inversion defined in Eq|2^. That is, 
exp(i2(/)r), Q"** transform into themselves under the Z2 
transformation and carry zero Z2 charges. 

In the next few sections we will apply the mapping we 
obtained so far to explore some aspects of the BEG of 
spin one atoms with antiferromagnetic interactions. We 
start with the symmetry broken states. 



The symmetry broken states of the BEC support very 
fascinating topological excitations, because of the pecu- 
liar structure of the internal space. In the next subsec- 
tions, I will explore these excitations in some detail, em- 
phasising on the influence of the 5^ and Z2 symmetry in 
the spinor BEC. 



A. Topological defects in BECs 



The Lagrangian admits a symmetry broken solution 
when the two-body scattering is weak. Topological de- 
fects are of particular interest because of the symmetry 
in Eq.|l^. They also play very important roles in the 
quantum symmetry restoring and are vital for the com- 
plete understanding of thermal and quantum phase tran- 
sitions. I will restrict myself to the static defects which 
can be directly detected in an experiment. 

In the previous section, two parameters, exp(ix) and n 
are introduced for the description of BECs. The ground 
state is degenerate under a gauge transformation of the 
U(X) group or a rotation of n on a unit sphere S'^. Fur- 
thermore, following EqjlJ, the wave function is indistin- 
guishable under a 180" n rotation and a gauge transfor- 
mation: X ~* X + ^- Therefore, the internal space for 
the order parameter is 7?. = [S*^ x S'^]/Z2; 5"^ x S"^ fac- 
torized with respect to Z2 is the product space of a unit 
circle and a unit sphere with two diametrically opposite 
points identified. That is (Xj n) = (x-I-tt, — n). This 
state should be defined as a quantum spin nematic state 
(QSNS). 

The symmetry found in Eq.|l^ was not recognized in 
previous works and the internal order parameter space 
was misidcntified as U{\) x S"^ in [Q, instead of [C/(l) x 
S'^]/Z2 discussed here. On the other hand, for a classi- 
cal nematic liquid crystal, the order parameter space is 
a unit sphere factorized with respect to Z2. The new 
complication in the QSNS, which originates from the ex- 
tra unit circle characterizing the superfiuidity, plays 
a very important role in topological defects. 
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A' X 




A' X 



It 27t 




(b) 





It 2jt 



(c) 




A' X 



It 2it 




(d) 




A' X 



It 271 




(e) 



Fig.3. Examples of paths in TZdS^ x S'^]/^2) and their 
projections in and 5*^. [S^ x S'^]/Z2 is represented 
by a torus or a cylinder with ends identified. The axis 
of the cylinder stands for a unit circle, along which x 
increases from to 2tt. Each disk, perpendicular to 
the axis, with boundary identified, is homotopic to a 
sphere; the boundary and the center of the disk are 
identified as the northern and southern pole of a two- 
sphere. Points such as A and A', which are diametri- 
cally opposite in the product space, are identical. The 
examples of paths homotopic to nonzero elements of 
X S'^]/Z2) are shown in a) (m = 1, n = 0) (see 
b) (m = — 1). A path homotopic to zero 
element is shown in c) (rri = 2, n = 0). Paths in d) and 
e) are homotopic to paths in a) and b) respectively. 

The linear or point defects in this case are determined 
by the classes of equivalent maps from a circle or a 
sphere of a real space into an internal space TZ, follow- 
ing the general principle for the classification of defects 
in a symmetry broken state | |58| , p9| , |60t . The ensemble 
of these classes, or the fundamental and second homo- 
topy group of space TZ are 7ri,2(7?-). Finally, textures. 



which correspond to defects with no singularities, are 
determined by a mapping from a three-sphere, a three 
dimensional space with surface identified, to TZ. Cor- 
respondingly, the spin configuration is homogeneous at 
infinity. The homotopic group for textures is n^lTZ). 



1. Z-i strings 



The fundamental homotopy group of space [S^ x 
S'^]/Z2 is a product of an integer group Z and a two 
element group Z2- The elements of a Z2 group are the 
indices for spin disclinations and Z is the winding number 
of superfiuid vortices. A closed path F belonging to the 
nontrivial element of the Z2 group is the one connecting 
two diametrically opposite points (%, n) and (x -I- tt, — n) 
in the space [S^ x S"^]. All other paths can either be de- 
formed continuously into this one or into topologically 
trivial ones of Z2 through escaping in a third dimen- 
sion(see Fig. 3). 

The corresponding wave function of linear singularities 
{Z2 strings) can be written as 

= n,$({?a}) 



01 



01 



2 Kq - 4o| 



m/2. 



hm n(e) - i?e(- 

w{m)+n 
hm Vclt) — — ; 



le-Col 



(27) 



Here ^ x + iy; m,n are integers and w{m) = 1/2 for 
odd ni and w{m) = for even m. We assume lines lo- 
cated at ^ = ^0 = xa+iyo. Each string is characterized in 
terms of (rn,n). is defined as the superfiuid velocity. 

The u;(m) dependent part of is present to ensure 
the single valuedness of the condensate wave function 
under the spin rotation. The spin wave function changes 
its sign when an atom moves along a Z2 string given in 
Eq.^, following the identity u{n)v{n) — ~u{—n)v{~n). 
The final spin wave function differs from the initial one 
by a minus sign. In a string (1,0), Vs of a half vortex is 
present solely to compensate the tt phase under n ^ — n 
rotationally in a Z2 string, as indicated in Fig. 4. A nc- 
matic Z2 string (l,rt) therefore corresponds to a tt spin 
disclination superimposed by a vortex with n + 1/2 fiux 
quantum, as shown in Fig. 4 p2[ |. 

This composite structure of linear defects is uniquely 
associated with the coherence of the condensate and is 
absent in a classical nematic liquid crystal. From an 
energetic point of view, a bare 7r-disclination carries a 
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cut along which the phase changes abruptly from tt to 
27r. This cut starting at the disclination ends only at 
the boundary of the coherent BEC and costs an energy 
linear in terms of the size of the system. For a similar 
reason, the energy cost to separate a half vortex and a 
TT-disclination at distance L is linearly proportional to L. 
The composite strings shown in Fig. 4 should be consid- 
ered as a result of the confinement of 7r-disclinations and 
half vortices in the spinor BEC. 




(a) 




Fig. 4 Z2 strings: a) spin wave functions along a loop 
enclosing a defect (see also the caption of Fig.l); b) tt 
spin disclinations (defined in Eq.^) superimposed with 
superfiuid vortices with half flux quanta (circular lines). 
In b), solid lines with arrows represent variation of n in 
defects. 

The energy of a linear singularity m — 1 , 71 ~ is 
pL / AM \n{L / a) , with L as the system size. We should 
emphasis that the spin component of the wave function 
carries a zero supercurrent. And is only determined 
by the polynomial of degree n + w{'m) in Eq.p7[ 

(— 1, n), (±3, n), (±5, n), ... linear singularities can be 
obtained from (— l,n) by continuous mappings(see ex- 



amples in Fig. 3). For instance, (— l,n) can be deformed 
into (1, ri) because of "an escape in a third dimension". 

For an even number m, n follows a closed loop on a 
unit sphere as a spin moves along the linear defects. The 
Berry's phase caused by the spin rotation along the linear 
defect is zero in this case, following an identity 



Im i dn ■ 



dfl d 
Att dn 



d 

u*{n)v*{n) X — — u(n)i;(n) = 



dn 



(28) 



The integral is carried over the area S on the unit sphere 
of n bounded by a path along which n varies as one moves 
along the defect. 

So there is no superflow of a half integer circulation su- 
perimposed when m is an even number. And vanishes 
in defects (m, 0) with an even integer m. More important, 
m = ±2, ±4, ... can be deformed into m = configura- 
tion as shown in Fig. 3 and are homotopically identical. 
There is only one family homotopically distinct nematic 
spin disclination (l,n). 

It is particularly interesting to have a closed loop 
string. The far field wave function of a closed loop linear 
defect can be expressed in a compact form as 



^({^4) = n„<i>({ej)u(n({e„ - ^o})Hna^c. - eo})), 

(29) 



lim n{^,9) = Re{-^)^^^^e„ 



hm v,(e,0) 



Here ep , are unit vectors in cylindrical coordinates 
(z, p, 6)] and ^ = p + iz. The circle S, — S,o is the center 
of the linear defect. The spin rotation in a constant 9 
plane leads to a half vortex; for a closed loop in Eq. |2^, 
the half vortex line forms a closed ring along eg direction 
and has no effect in the far field. 

The stability of a closed loop is a rather subtle issue 
and deserves some attentions. It depends on competi- 
tions between the energy associated with the tension of 
the loop and the spin fluctuations which unfavor a singu- 
lar structure. Under certain conditions, one indeed can 
show that a close loop structure is stablized |6|] . 




(a) 



(b) 



Fig. 5 The far flelds of a closed loop 7r-disclination are 
identical to those of a hedgehog. The topological charge 
of Pontryagin type spreads over the entire loop while in 
a hedgehog, the charge is localized at the center. 
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However, if the loop shrinks, one ends up with a 
monopole which I am turning to in the next subsection. 
The intriguing connection between a closed loop tt discli- 
nation and a hedgehog also lies behind a recent inter- 
esting work [Q. In the non-Abelian field theory anal- 
ogy of the situation discussed here, it was pointed out 
that a monopole could exist in a form of an Alice string. 
One can examine the connection field (defined in section 
VI) of a closed loop to show indeed the Pontryagin field 
strength is nonsingular at the center while the overall 
field strength at a large distance is still of a unit charge. 
The small distance behavior is distinct from a point-like 
defect. It implies that the unit charge of the Pontrya- 
gin type is uniformly distributed along the tt disclina- 
tion. In a large loop limit, the linear charge density is 
vanishly small but with an integrated charge of unity. A 
charge exists as a closed loop string, with no singularities. 
Thus, the BECs of "^^Na might also be playgrounds for 
the physics of Alice strings. A close analogy of the Alice 
string physics was found in |^5| where a hypermonopole 
is transformed into an antihypermonopole when moved 
along a 7r-disclination. 

2. Hedgehogs 

The second homotopy group Tr2{[S^ x S'^]/Z2) has non- 
trivial elements. Point defects are of the form of "hedge- 
hogs" with the degree of mapping given by N^, 

In a cylindrical coordinate (z,p,6), the far field wave 
function of a Nh = 1 hedgehog located at the origin is 

ii?e|^ cos(0„)(i;2 - u^) + ^i?e|^ sin(0,)(«2 -|- u^)] 

lim n(^, e) = i?e— Gp + Im—e^, 

lim v,(e,0) =0. (31) 

where = p + iz. 

However, consider a hedgehog with an index Nh- One 
can always move n(r) along a closed path F in space 
[S^ X S'^]/Z2, which is homotopic to the nonzero ele- 
ment of 7ri([S'^ X S'^]/Z2)- This connects (x(r), n(r)) into 
(x(r) -I- TT, — n(r)) and transforms a hedehog with index 
Nh into a —N^ hedgehog. In general, under the influence 
of 7ri([S'i X S'2]/Z2), the elements in tt2{[S^ x S'^]/Z2) can 
be transformed into one another (from A^^ into —Nh) by 
moving (x, n) along the path F. This was emphasised 
in the investigation of a classical nematic liquid crystal 
in where TZ = Physically it can be done by 



introducing a tt spin disclination superimposed by a tt 
vortex and moving the hedgehog along the linear defect; 
thus a linear singularity given in EqJI^ defines a contin- 
uous deformation from Nh into —Nh- The point defect 
is characterized only by the modulus of A^^i and a hedge- 
hog with an index A''^ is homotopically identical to that 
with —Nh- Following this discussion, the energy barrier 
involved in this deformation is the energy of the linear 
defect and is linear in terms of the sample size. 

Despite the homotopical indistinguishability in the 
presence of linear defects, the positive and negative 
hedgehogs can be physically distinguished because of the 
coherence of the BEG, unlike the situation in a classical 
nematic liquid crystal. If a tt disclination is introduced 
alone, by moving a hedgehog around the linear defect, 
the condensate will acquire a 7r-phase with respect to the 
original one. The 7r-phase difference, which can manifest 
itself in a Josephson type of effect, is one of the signatures 
left behind by the positive or negative hedgehogs. 

The second way to distinguish the positive and neg- 
ative hedgehogs is to look at the local connection field. 
The idea here is to introduce a spin-^ collective excita- 
tion, which carries a half charge with respect to the con- 
nection fields. Unlike atoms which carry zero charges, 
in a positive hedgehog configuration, the spin-^ object 
experiences a connection field of an opposite sign com- 
pared with that of a negative hedgehog. This results in 
the Berry's phases of different signs in the positive and 
negative hedgehogs. 

The other influence of the Z2 symmetry on the point 
defect is that a hedgehog is identical to a closed loop Z2 
string as illustrated before. Though the spin hedgehog 
and the closed loop tt disclination are topologically iden- 
tical, when spin fluctuations are insigniflcant the closed 
loop string has a higher energy because of the linear ten- 
sion of a string. In this limit, we expect a closed loop, 
once created, collapses and leaves a spin hedgehog be- 
hind. The energy of a spin hedgehog is proportional to 
pL/AM, linear in terms of the system size L. Some static 
and dynamical aspects of the hedgehogs have also been 
investigated recently in |Q . 

The homotopical identity between a 7r-disclination and 
a hedgehog, and the conversion between the positive and 
negative hedgehog in the presence of 7r-disclinations are 
fascinating properties of the quantum spin nematic BEG. 
The Z2 symmetry, after all, influences every aspect of the 
point defects under the consideration. 

At this moment, we want to point out that the "hedge- 
hog" in a ferromagnetic BEG is a " spin hedgehog" super- 
imposed with a superfiow. Along ~z axis, the circulation 
along eg is quantized at / dl • Vx = 1, representing a vor- 
tex line ending at the monopole or a vorton discussed in 
| ]65| , p9[ . These features are absent in the QSNS for the 
reason that the Berry's phase as shown in Eq.BS is zero. 
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3. Hopf textures 



B. The effects of an external magnetic field 



Textures are characterized by a Hopf mapping from 
to [5*^ X S'^]/Z2 in the presence of the supcrfluid compo- 
nent. In this case 7r3([S'^ x S'^]/Z2) = Z is an integer 
group. The Hopf number which characterizes the degree 
of mapping can be written in terms of n as 



In the presence of an external magnetic field along z 
direction, r2i.2 = nzt e^gt-iBH/Ac-zp, following Eq.pd]; 



H = — d-'rA • V X A, 
Sir J 

V X A = e'''"'eijknadjnbdkncei. 



(32) 



In a cylindrical coordinate (p, z,9), the p — z plane is 
homotopical to a two-sphere, with the northern pole iden- 
tified as the center of the texture and the southern pole as 
the boundary. The configuration on a p— z plane is there- 
fore homotopically identical to a mapping from 5*^ to a 
[S^ X S'^\IZi and can be considered as a Skyrmion on the 
p—z plane. By rotating it along z axis, we obtain the tex- 
ture with Hopf number H = 1. The center of the texture 
is located along a circle = around z axis = p + iz). 
Similarly, surface /m(f — ^o)/Rg{^ — ^o) — const with 
boundaries identified is also homotopically identical to a 
two sphere S'^. Under the Hopf mapping, the spin con- 
figurations on these planes are the same as Skyrmions. 

So in a cylindrical coordinate (z, p, 6), the far field wave 
function for a texture with H = 1 is 

^aej) = n„<i)(Ua}) 

[^sing(|e„ -^ol)i?e(J^i^exp(^0„))(^;2 - w2) 

+ ^ sin.g(|ea - eo|)/m( exp(i0„))(z;2 + u^) 

+ COS g{\^a - £,o\)uv] 

hm n(e,0) =sing(|e-Co|)[i?e|^ep-H/m|:i|^e,] 

+ cosg{\£,- Co\)e,, 
lim v,(C,0) = 

^ — >oo 



(33) 



Here g{\£_\ = 0) = and g{\£^\ = oo) = tt. Once again, 
the superflow is zero in a Hopf texture in the QSNS 
by contrast to that in a ferromagnetic BEG. 

Since the order parameter is homogeneous at the 
boundary of the space, the texture can considered as a 
"particle" of a finite extend <^o- The energy is estimated 
as p£,o/AM and decreases when the texture shrinks; this 
implies an instability towards collapsing. It can be en- 
ergetically stablized only when high derivative terms are 
taken into account, as suggested in [Q. However, this is 
beyond the validity of the nonlinear sigma model where 
only the long wave length dynamics is correctly charac- 
terized. 



gp.BH 
1 = — ; Bz 

4C2P 



(34) 



The Lagrangian is the same as that without an external 
field (up to a constant term). However, the constraint 
that n has to satisfy depends on the external field. In 
the absence of an external magnetic field, Eq.9 is auto- 
matically satisfied for 1 = n x 9tn/4c2. Cs is of the form 
of the 0(3) NLaM. An external field breaks the 5'^ sym- 
metry and confines the low frequency sector of n in a 
plane perpendicular to H itself, i.e. 



n • = 0. 

4C2P 



(35) 



The Lagrangian in the presence of a magnetic field is 
that of an 0(2) NLctM; it has the 5^-symmetry at the 
frequency to pbH . At the high frequency limit, fol- 
lowing the equation of motion, n precesses in a field 
4c2l(a;), much larger than the external one H and the 
5^-symmetry is restored. 

As a consequence, if an external magnetic field H — 
Hgz is applied, ng is in the xy plane and the parameter 
space for ng is a unit circle. The ground state is degen- 
erate under the group \U{1) x S^]/Z2, following Eq.l3. 
The order parameter space for the quantum spin nematic 
state therefore is 7?, = [S*^ x S^]/Z2, which represents the 
product space of two unit circles but again with two di- 
ametrically opposite points {x, n) and (x -I- tt, — n) being 
identical. 

The fundmental group 7ri([S'^ x S^]/Z2) has nontriv- 
ial elements. The corresponding wave function for linear 
defects can be written in the same form as that in Eq,27, 
TO, n are integers. However in this case, all defects (to, n) 
are topological distinct (at energy scales lower than the 
external field), m = ±1 and n = ± are positive and 
negative tt spin disclinations, shown in Fig. 2. 

So far we discuss the static properties of the bulk de- 
fects. The dynamical aspects of these defects as well as 
surface defects will be addressed elsewhere. We also no- 
tice that in the presence of the external field, by contrast 
to the classical nematic liquid crystals, n in the QSNS is 
not parallel to the external field. This results in sound- 
like collective modes even in the presence of an external 
field, but with number of modes reduced from two(at zero 
field) to one. 

We have restricted ourselves to the weak magnetic field 
limit and neglected the possible quadratic Zeemann shift 
Hqz = Yl,a QH^^za (t^^e external field H is along 
direction). An inclusion of the quadratic Zeemann shift 
yields an additional term Cqz = / pQH^{n^ + ny)dx to 
the NLcrM derived in section II. The main effect of this 



13 



contribution is to align n along the external field. There- 
fore, when this shift dominates, the ground state is left 
with a double degeneracy: n = and n — — e^. The 
spin wave develops an energy gap of the order QH^ . We 
will focus on the zero magnetic field case in the rest of 
discussions. 

(2jt,0) (2jt, In) 



(It, II) 





(0 0) (0,2it) 



(2ir,0) (2ir, 2it) 



(It, n) 
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(0 0) (0,2ii) 
(2it,0) (2it, 2it) 
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(0 0) (0,2it) 
(2it,0) (2it, 2it) 



(It, It) 



(0 0) (0,2it) 




A" X A, 




Fig. 6. Examples of closed paths in TZ — [S^ x S^]/Z2 and 
projections in and . Here the surface of a torus-2, 
or a square with boundaries parallel to each other iden- 
tified, stands for [S^ x S^]/Z2- The points such as (tt, tt) 
and (0,0), or (0,27r), or (27r,27r), or (27r,0) are identi- 
cal, representing two diametrically opposite points in the 
product space S*^ x 5^. The examples of paths which are 
homotopic to nonzero elements of 7ri([S'^ x S'^]/Z2) are 
shown in a). (m=l, n=0), b). (m=l, n=-l), c). (m=-l, 
n=l), d). (m=-l, n=0). All are topologically distinct. 



IV. SPIN DISORDERED QUANTUM 
CONDENSATES 



For the quantum nematic state considered in the previ- 
ous sections, the full Hamiltonian includes Tis for the spin 
sector (S*^), Tic for the superfluid sector (S"^), the cou- 
pling between two sectors Tisc and the constraint that the 



wave function remains invariant under a local Z2 trans- 
formation, i.e., n(x), X + tt) = 5'(n(x), x). To study 
quantum nematic fluctuations or quantum spin nematic 
order-disordered transitions, one has to integrate out the 
phase degree of freedom (superfluid component) under 
the constraint and obtain a renormalized spin dynamics. 

While a full treatment of the quantum spin nematic 
fluctuations taking into account the influence of the Z2 
symmetry will be present in the next section, here I 
will be interested in a spin disordered condensate due to 
strong interactions between spin wave excitations in the 
0(3) NLctM; in this limit, the coupling between the phase 
fluctuations and the spin fluctuations is irrelevant. The 
procedure is carried out self-consistently: first we obtain 
results based on the o(3) NLcrM and then consider the 
influence of the sector and the Z2 symmetry. The re- 
duction from [5^ X S'^]/Z2 to 5^ for the spin dynamics is 
possible based on the following considerations. 

First, the phase fluctuations in the BEC and the ne- 
matic order fluctuations can affect each other via a di- 
rect coupling Csc- Nevertheless, the inclusion of Use 
only results in higher derivative terms; the coupling 
thus becomes important only at a rather high frequency 
hp^/'^ /2M, which is also the upper cut-off energy of the 
NLcrM. 

Second, the constraint in Eq.^ leads to another cou- 
pling of a pure topological origin; for instance in the static 
limit, it leads to superimposed Z2 strings as discussed in 
section III. The spin nematic order-disorder transition 
depends on the Z2 gauge fields and spin waves in the 
S"^ sector. However, if the phase fluctuations are weak 
across the nematic order-disordered transitions driven by 
the zero point rotations of each atom, the Z2 fields are 
effectively frozen and irrelevant to the discussion. 

Indeed, in 1 -I- 1 dimension discussed in the next sec- 
tion, the nematically disordered state exists in the weakly 
interacting limit where the phase fluctuations are still 
negligible. In terms of quantum tunneling, the nematic 
disordered state in 1 -I- 1 dimension is driven by tex- 
ture instantons. Quantum tunneling of the Z2 strings 
(or Z2 instantons) turns out to play little role in the ne- 
matic order-disordered transitions because of the coher- 
ent phase sector. This again implies the Z2 gauge flelds 
are not relevant to our discussion. In 3 -I- 1 dimension, 
this is also the case and the transitions are not driven 
by the Z2 gauge fields if the phases of the condensates 
remain coherent around the nematic order-disorder tran- 
sition point. 

As in other symmetry broken states, quantum fluctu- 
ations in n exist in the spinor BEC because n doesn't 
commute with the Hamiltonian and is not a conserved 
quantity. The quantum fluctuations in the QSNS can be 
studied by considering the NLcrM, 



(36) 
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Here 



A. Rotation symmetry restored correlated states in 

3d 



/ = (167r)i/2(pAa3) 1/6, ^ (37) 

in a single trap limit. We also introduce dimensionless 
length and time: f = rp^/^^ f = rvsp^^'^^ and Vg — 
^4c2p/M. Derivatives are defined as {df,dx,dy,di). 

In the typical low frequency fluctuations or spin waves, 
|L| is much less than unity and each individual atom 
remains to be in a spin state M(n(r))w(n(r)), but with 
n(r, t) varying in the space and time. The quantum fluc- 
tuations of the nematic order can be estimated in the low- 
est order approximation in the weakly interacting limit. 
This is similar to the spin wave approximation made for 
Heisenberg antiferromagnetic spin systems (HAFS) pTj . 
The correlation function of 6Q°''^ = Q°'^- < Q"^> 
(a = X, y, z) in the QSNS can be calculated as 

<5g"'3(r,i)5Q"''=''(0,0)> = 
n 

+SoaSof3'Pf3rjPa'7i + (^Q/Jf^Oa' Pq,,P/3',, )C,, (r, t) . (38) 

<> represents an average over all configurations. =< 
(5n^(r, t)(5n^(0, 0) > reflects the dynamics of the "direc- 
tor" ; T] — ±1 label eigen modes and labels the direction 
of no = Gz- Pari = Gq, • 6^, whcrc is the unit vector 
of the rjth eigenmode; for ng = e^, = (e^; +iriey)/^/2. 
We assume (5n = n — no is small and different modes 
interact weakly. Following Eqjs^, 

/ N „ f ,^ , expfipi — iiut) 
C^r, t) = f dPdu;-^ 39 

cJjP are dimensionless variables. The amplitude of C,, is 
proportional to the parameter /. Eq. ^ is valid when 
/ is much less than unity. When / increases, spin wave 
excitations start to interact strongly and Eq.^ becomes 
invalid. 

/ is a measure of the amplitude of quantum fluctua- 
tions in the BEC of '^^Na. Following EqjH the energy of 
the system consists of two parts, a) The potential energy 
?i^p(Vn)^/2A/, which is the energy cost in the presence 
of a slow variation of n, aligns n of different atoms. It de- 
termines the bare spin stiffness. And b) the zero point ki- 
netic(rotation) energy ?i^P/2/o, where Iq = {4:C2p/h'^)^^ 
can be considered as the effective inertial of an individ- 
ual atom. This rotation energy originates from the two- 
body scattering in the microscopic Hamiltonian and the 
inertial is inversely proportional to the scattering length. 
The zero point rotations tend to disrupt the order of n 
between different atoms. is a square root of the ra- 
tio between the potential energy at an interatomic scale 
?i^p^/^/2m and the zero point kinetic(rotation) energy 
h?' /2Iq of an individual atom. 



At zero temperature the o(3) NLctM in Eqjs^ has or- 
dered and disordered quantum phases at d > 1, depend- 
ing on the parameter /. The renormalization group (RG) 
equation ffow is determined by the interactions between 
collective modes. The scatterings between these spin 
waves are determined by the local nonlinear spin dy- 
namics, especially spin waves around some slowly varying 
spin configurations. In the following, I will take this RG 
equation point of view of the NLtrM: / is the quantity 
determining spin correlations in the ground state of the 
BECs. 

The RG equation for / can be obtained by first inte- 
grating out fluctuations within e~' < IpI < 1, e~' < e < 1 
and then rescaling p — > pe' , e — > ee' . In 3 -I- 1 dimension, 
the RG equation takes a form 

f = -2/(/c-/) (40) 

Within the framework of the NLctM, fc = Stt^ in d = 
3(see Appendix B). However, the exact value of fc is de- 
termined by the details of the short-range behavior of the 
system and should be obtained only by numerics. The 
state where the long-range nematic order is absent due to 
zero point motions of collective variables is called a " spin 
disordered quantum condensate" (SDQC), to be distin- 
guished from the "QSNS" where the long-range order is 
present. 

At the strongly interacting or a high density limit 
/ > /c, the zero point kinetic energy dominates and 
the spin stiffness is renornialized to zero at a long wave 
length limit. Especially, in an extremely quantum disor- 
dered phase where the scatterings between two atoms are 
strong and the inertial is small, n of each atom rotates 
independently and the spin wave function of atoms only 
correlates at an interatomic distance, n of each spin ex- 
periences fast diffusive motion on the unit sphere S"^ with 
a mean free time l/2c2p much shorter than 2M /h? p^/^ , 
due to the strong zero point rotation. The excitation 
spectrum has a gap of order of 2c2p. 

As the two-body scatterings get weaker, the inertial of 
each individual atom /o becomes higher and the motion 
of n on the unit sphere is slower because of a lower zero 
point rotation energy. Thus, n of different atoms starts 
to correlate at a flnite length ^, which is a function of the 
parameter /, much longer than the interatomic distance. 
As this happens, the resultant effective inertial increases 
proportional to the number of atoms within the range 
of the correlation length, /(^) ^ ^'^/2c2. And the spins 
of different atoms precess a collective temporal motion 
on the unit sphere with the mean free time r. Follow- 
ing Eq.||, L(^) ~ /c2T and the total rotation energy 
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is L(^) /2/(^). This has to be balanced by the poten- current experimental conditions. For pAa^ 



10- 



as IS 



tial energy p£,/M, leading to r ~ ^/ y^C2p/M . It implies 
that in this spin correlated state, collective excitations 
propagate with a bare spin wave velocity ^CipjM up 
to a scale ^. Finite-range nematic order exists in this 
case and the correlation length can be estimated using 
the renormalization approach. In the current situation, 
^ = p^^/'^/c/[(pAa^)"'^/^ — /c]. Local nematic order is es- 
tablished up to a distance ^ in the BECs but the state is 
rotationally invariant. A phase transition takes place at 
pAa^ = 8 X 10"^ X /6. 

At weakly interacting or a low density limit f < fc, 
the potential energy dominates and the spin stiffness 
flows to a finite value under the renormalization group 
transformation. Long-range order occurs. All spin one 
Bosons rotate as a rigid body with an effective inertial 
I{N) — N/c2P- The energy gap in the excitation spec- 
trum decreases accordingly 



gap 



2I{N) 



h'^C2P 

2N 



(41) 



as the number of Bosons N increases and vanishes as 
a thermal dynamic limit is approached. The rotation 
symmetry is broken in the ground state of the BEC as 
nematic long-range order is established. 

We want to make the following remarks concerning the 
statements made so far; 

1) . The quantum spin order-disorder transition in 3d 
is driven by strongly interacting short wave length ex- 
citations. Therefore, in principal, it is no longer valid 
to neglect the coupling between spin and superfluidity as 
we did for the low frequency sector. We, however, believe 
the picture present here is a qualitative correct one and 
the coupling can at most shift the value of fc- 

2) . In the weakly disordered limit, the linear tension 
of a Z2 string at a scale L is of order Lp~^/'^Ecore, with 
the core energy Ecore of order p^^^ /2m. Large loop Z2 
instantons are extremely rare and can be neglected. This 
justifies the reduction we made at the beginning: The Z2 
gauge fields fluctuate weakly in this limit. 

3) At a high density limit, one should also take into 
three-body, four-body scatterings. This can further mod- 
ify the short distance dynamics but will not affect the 
conclusions arrived above in a qualitative way. 

4) One should be cautious about the definition of 
"phase" since the alkali atoms under investigation are in 
a long lived metal stable gaseous state. The life time of 
atomic gas is limited by three-body inelastic collisions . 
The collision rate increases dramatically as the density 
increases. For sodium atoms, the life time of the metal 
stable gas is given as Q] = 10~^^ x (pAa^)~^ x lOOsecs. 
If a reasonable experiment measurement can be per- 
formed at a time scale 10ms, then p/S.a^ should be less 
than 10^**. This sets a limit on fc in order for the quan- 
tum disordered nematic phase to be probed under the 



(42) 



in the experiments, long-range nematic order should be 
observed. 



B. Rotation symmetry retored states in Id 

Quantum fluctuations of collective variables are most 
prominent in BECs confined in highly anisotropic traps. 
The N spin one bosons ground state wave function in 
this limit does not live in the zero momentum sector of 
the Hilbert space. The NLcrM turns out to be the most 
powerful approach to capture spin correlations in a ro- 
tationally invariant ground state. For simplicity, we ap- 
proximate highly anisotropic traps as Id nematic BECs. 

In Id traps at zero temperature and zero field, one 
should expect there will be no long-range order and the 
state is nematically disordered, following the renormal- 
ization group results of the NLctM |^ . The RG equation 
in this case can be shown as 

dl 2tt' 

which always flows into a strong coupling fixed point (dis- 
ordered state) at the low energy and long wave length 
limit. The correlation length is ^ = exp(47r//) [p is 
the linear density). That is, in a Id trap, the excitation 
spectrum should have an energy gap following ]4l| ]. 

In the Id spin nematic BEC under consideration, 
where the long wave length spin wave fluctuations are 
most significant and 1(2;) 0, Tisc can be effectively 
taken to be zero. To a good approximation, the phase 
(or density) fiuctuations and the nematic spin fluctua- 
tions are indeed decoupled (the Z2 gauge fields be also 
frozen as discussed at the end of this section). Therefore, 
the ground state can be spin disordered but phase coher- 
ent. The Id BECs of "^^Na have long-range order in x 
but short-range correlations in n in the weakly interact- 
ing limit. 

These spin disordered states mimic the quantum spin 
liquid states proposed in the literature of Heisenberg an- 
tiferromagnetic systems (HAFS). They represent nematic 
liquids, where n of "^^Na atoms in the BEC aligns toward 
a certain direction up to a length scale ^ and the order 
is disrupted afterwards due to the zero point motions 
originating from the spin-dependent two-body scatter- 
ings. The zero point motions in n also lead to quantum 
tunneling of spin textures, with the tunneling rate per ^ 
inversely proportional to the energy gap. 

We notice that the NLctM obtained from the micro- 
scopic Hamiltonian in this paper doesn't have a 0-term 



La = 



dn dn 

drdxn • ^ x — , 
or ox 



(43) 



which is generally present in the HAFS studied before 
pl[ . The presence of a 6'-term in the NLcrM can dramat- 
ically affect the RG flow. A 9-teim with 9 = tt, as in 
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the half- integer HAFS, results in an additional term in 
the RG equation In fact, the additional term aris- 
ing from a 0-term forces the RG to flow into a massless 
fixed point in Id and implies a power law decay of the 
correlation function or gapless excitations. 

The absence of a 6'-term in the NLtrM derived here, 
which prevents the RG from flowing into a massless fixed 
point and ensures an energy gap in the excitation spec- 
trum of the ncmatic state, can also be justified on a 
ground of the Z2 symmetry. It follows a general fact that 
the Berry phase in an ordered state vanishes as a result 
of cancellation between u{n) and v{n) contributions, as 
shown in Eq.p8|. 

We should note that in a nematic state, the textures 
with 1 and —1 indices in 1 + 1 dimension are homotopi- 
cally indistinguishable under the influence of the Z2 sym- 
metry. This is similar to what happens to the point de- 
fects in 3d: the spin hedgehogs with index Nh and ~Nh 
are indistinguishable as emphasised in section IV. 

The Z2 symmetry in our theory also leads to tt spin 
disclinations superimposed with half vortices in 1 + 1 di- 
mension: Z2 instantons. Unlike texture instantons, the 
action for an individual Z2 instanton is logarithmic diver- 
gent. When the textures are already present, the interac- 
tion between the spin disclinations in the Z2 instantons 
becomes finite ranged due to the loss of the spin stiffness 
at a distance longer than ^. However, the half vortices in 
the Z2 instantons still interact with each other via a log- 
arithmic potential when the phase remains coherent. So 
the quantum tunneling of the Z2 instantons are absent in 
disordered nematic states in the weakly interacting limit. 

The suppression of the Z2 instantons caused by the 
half vortices also suggests that the Z2 gauge fields, when 
coupled with a phase coherent BEG, be "frozen" in 1 -I- 1 
dimension in the weakly interacting limit. The Z2 gauge 
field is effectively decoupled from the spin dynamics and 
the spin order-disorder transitions are driven by the spin 
waves in the sector, which is self-consistent with the 
assumption we made at the beginning of the section. On 
the other hand, when both the phase and the spin sectors 
are disordered, we expect the presence of the Z2 instan- 
tons in 1 -|- 1 dimension. A quantitative calculation of the 
multi-influence between the it spin disclinations and the 
textures, especially the tt spin disclinations effect on the 
Haldane gap is still absent fldj . 

V. SDQCS UNDER THE INFLUENCE OF Z2 
FIELDS 

The influence of the Z2 symmetries in symmetry un- 
broken states, say SDQGs, can be most conveniently ex- 
amined in the lattice gauge field model discussed in sec- 
tion 11. The action in Eq.|2^ represents an model and 
a Heisenberg model coupled with a Z2 gauge field cr^r'- 
In the limit when K goes to infinity, the Z2 gauge fields 



are frozen and cr^r' = 1; Eqj2^ reduces to an model 
for the phase sector and a Heisenberg model for the spin 
sector. The spin dynamics in this case is identical to that 
of an NLcrM, as discussed in details in an early section. 
The dimensionless coupling constants for the XY model 
and the NLcrM in optical lattices are 

/;■'• = /f'--/^ (44) 

where g and u as defined in Eq.|l^ depend on the spin 
susceptibility and the compressibility. By carrying out 
calculations similar to those in section 11, one obtains 
u = 47rp(2a2 -I- ao)/3MNo and g = 47rp(a2 - ao)/3MNQ 
in terms of scattering lengths 09,27 the number density p 
and the mass of atoms M. 

In the weakly interacting limit where /f '' is much 
less than unity, the rotation symmetry is broken in the 
ground state and nematic order occurs. This order leads 
to birefringence of a polarized light and enhanced 
small angle light scatterings. However, when /° '- is 
much larger than unity, the interaction between Gold- 
stone modes results in a spin disordered liquid. Similar 
to a resonant-valence- bond (RVB) liquid in Heisenberg 
antiferromagnetic systems, a spin disordered liquid is ro- 
tationally invariant and thus is free from birefringence. 

When J'^-'^ approaches zero, Eq.p3| becomes identical to 
a pure Z2 gauge theory, which is dual to an Ising model 
II72I . There are two phases for a pure Z2 gauge theory at 
the spatial dimension greater than one. When K > Kc, 
the Wilson loop integral (defined in the next section) de- 
cays exponentially as a function of an area enclosed by 
the Wilson loop, i.e. in an area law. When K < Kc, the 
Wilson loop integral decays exponentially as the func- 
tion of the perimeter of the loop. In the former case, two 
Z2 charges interact with a potential linearly proportional 
to the distance between charges, therefore are confined. 
This has a rather profound impact on the quantum num- 
bers of quasiparticles which will be studied in some detail 
in section Vll. 

In general, the BECs characterized by Eq.^ exhibit 
extremely rich physics because of a coupling between 
the spin degree of freedom and Z2 gauge fields. The 
influence of the Z2 symmetry in symmetry broken states 
was discussed in the previous section and here we want 
to further explore this influence in some symmetry par- 
tially restored states. Particularly, we will be interested 
in fractionalization of topological excitations due to the 
Z2 gauge fields. 



A. Fractionalization of integer vortices 

There are three types elementary topological defects in 
the BEG of 23 iVa ||: 
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1) a half-vortex around which the phase of the conden- 
sate changes tt; 

2) a TT spin discHnation around which n rotates by 180°; 

3) a frustrated plaquette with Yin '^ij — o^' ^ ^2 
vortex. A Z2 vortex emits a string of aij = — 1 hnks ter- 
minated either at the boundary or at another Z2 vortex. 

In general, when the U{1) symmetry or the S*^ symme- 
try is broken, an individual half vortex or a 7r-disclination 
or a Z2 vortex is an energetic catastrophe because of a 
cut, running from the defect to infinity, along which the 
phase of the condensate has to change abruptly. As a re- 
sult, the Z2 vortices have to be confined and form bound 
states with either the half vortices or the tt disclinations 
in the phases where the U{1) or and the Z2 symme- 
tries are broken. In a polar state studied before, where 
all three symmetries are broken, only composites of the 
half vortices and the tt spin disclinations superimposed 
with the Z2 gauge vortices are free objects |37| , ^ . On 
the other hand, in the states where both the U{1) and 
symmetry are unbroken, the Z2 vortices are free ex- 
citations. 

To demonstrate fractionalization of topological excita- 
tions in the SDQC explicitly, we consider a specific sit- 
uation where the U{1) symmetry is broken but the S'^ 
symmetry is restored. This happens in Id in a weakly 
interacting limit and also in higher dimensions. So the tt 
disclinations are always deconfined and are liberated. By 
integrating out n of the nematic spin liquid, one obtains 
an effective action for exp(ix). 



S 
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with JxY, being finite, renormalized by gapped spin 
fluctuations [Q. Eq.^ represents an XY model cou- 
pled with a Z2 gauge field. A gauge model similar to 
Eq.^ was proposed for the classical 0(3) nematic order- 
disorder phase transitions fz^ j . A similar model was also 
considered and investigated in the context of cuprate 
high-temperature superconductors. 

The BECs described by the reduced action in Eq.^ ad- 
mit one ordered and two disordered phases. At large 
limit, a nematic order-disorder phase transition takes 
place, with integer vortices condensed but half vortices 
gapped. The energy of half vortices per unit length is 
finite in the disordered phase. This state as discussed 
in section VII supports spin one excitations and is a 
fractionalized phase. At small limit, another phase 
transition occurs, but with half vortices condensed: the 
energy of a half-vortex per unit length vanishes in the 
disordered phase. It is a confining phase where only spin 
two excitations exist. Along the Jxy = axis, there is a 
transition between these two disordered phases, the frac- 
tionalized one and the confining phase, as a pure gauge 
field phase transition. 

When Jxy is large, the U{V) symmetry is broken. Un- 
like in one-component BECs where an elementary topo- 



logical excitation carries a unit circulation, an elementary 
excitation in the BEC with a hidden Z2 symmetry is a 
half vortex. To see this, let us start with the limit where 
is large but finite. The energy density of an isolated 
half vortex per unit length is 



h.v. 



pirh / L L^ 



(46) 



with the second term being the energy of a cut. (L is 
the size of the system in xy plane). Moreover, two half- 
vortices are bound to form an elementary excitation with 
a unit " flux quantum" . 
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Fig. 7 A Z2 vortex emits a string of 
resented by short-thick lines. 



-1 links, rep- 



However, when a Z2 vortex as shown in Fig. 7 is in- 
serted, all arr' along the cut have negative signs and the 
energy catastrophe of a half vortex is removed. The en- 
ergy of a Z2 vortex is from the core where Yin '^ij ^ 
negative sign there. The finite energy cost of the gauge 
vortex is proportional to i^^j (assuming Z2 is deconfin- 
ing). The total energy of the composite in the unit of 
pTrfi^/SM is just InL/o, up to a logarithmic acuracy and 
is one fourth of the energy of an integer vortex. There- 
fore, an integer vortex will fractionalize into two half vor- 
tices and a half vortex becomes an elementary object. Let 
us emphasis this arises because of the discrete Z2 gauge 
fields and doesn't happen in spinless BECs or an XY 
model. 



B. BECs of singlet pairs 

Upon taking into account other higher order processes, 
we have to include a hopping integral of singlet pairs as in 
Eq.p3|. The pair hopping integrals J^^'^^p g^j.g particularly 
important when g is large and states with odd number of 
atoms in individual wells are energetically costly: they 
are not allowed to have L = state of a rotor due to 
parity constraint and have a higher rotational energy. It 
is energetically favorable to have even numbers of atoms 
on each well arranged into singlet combinations, which 
can be interpreted as binding of atoms into singlet pairs 
with a pair binding energy 
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Ep{Q) = 2E{Q + 1) - E{Q + 2) - E{Q) 



(47) 



In the hmit when g — > cxd a crystal phase is possible 
for even numbers of atoms per well and corresponds to 
a spin singlet insulator. The superfluid phase is a con- 
densate of singlet pairs (SSC), in which tunneling of in- 
dividual atoms between the wells is suppressed and only 
singlet pairs are delocalized. The origin of pairing in 
this case is not the attraction between individual bosons, 
but a singlet formation on the scale of individual wells. 
This reminisces the "attraction from repulsion" mecha- 
nism of electron pairing proposed by Chakravarty and 
Kivelson for high Tc cuprates, Cgo, and polyacetylene 
iQ . The basic idea of this type is also the foundation of 
Anderson's interlayer tunneling mechanism for the high- 
Tc cuprates [75| ] . In general, when a liquid is one particle- 
incompressible but two-particle compressible due to ei- 
ther a spin gap as in our case, or more general due to in- 
frared catastrophes, the two-particle tunneling between 
liquids is encouraged and long-range order of pairs can 
be established (7^ . 

Formally speaking, a large pair-tunneling term reduces 
the C/(l) symmetry to the Z2 symmetry; only cj) — Q 01 tt 
on the unit circle is relevant for the low energy physics. 
The planar spins in Eq.^ therefore have to be replaced 
by Ising spins Sr ±1 and one obtains 



(48) 



The action in Eq.^ admits only two distinct phases, as 
pointed out by Fradkin and Shanker |7^]. Along the 
= 00 axis, a usual Ising order-disorder transition 
takes place. Along the Jj — axis, there will be a tran- 
sition from the Z2 confining phase to the Z2 deconfin- 
ing phase. But along the = axis, no transitions 
occur. So two symmetry broken but topologically dis- 
tinct phases are allowed: a) an Ising ordered phase with 

< S'i >^ and confining Z2 fields and b) an Ising disor- 
dered phase with < Si >= and deconfining Z2 fields. 

For BECs, the first phase corresponds to a conven- 
tional single atom condensate with < cosx and 
the latter one is an exotic singlet pair condensate where 

< cosx >= but < cos2x >7^ 0, first pointed out in 

Topological excitations have distinct composi- 
tion in these two symmetry broken states. The Ising or- 
dered phase breaks the Z2 and U{1) symmetries. The 
half vortices and the Z2 gauge vortices are confined. 
However, in the Ising disordered state, the Z2 symme- 
try is restored while the U{1) symmetry is still broken. 
The half-vortex and Z2 vortex are liberated and are free 
excitations. The transition between these two states is 
of the Ising universality. The other interesting aspect of 
the singlet pair condensate is the a.c. Josephson effect. 
For the pair condensates. 



X2) 



where xi,2 and ^1^2 are the phases and chemical poten- 
tials of two BECs. Note that the factor 2 appears in the 
right hand side of the equation for the pair BEC. A more 
realistic proposal was given in |Q . 



VI. HIDDEN TOPOLOGICAL ORDER 



As emphasised through out this paper, the spinor 
BECs of spin one atoms have many fascinating proper- 
ties associated with the unique topology of the spin order 
parameter and the entanglement of spin-phase degrees 
freedom. These interesting topological aspects result in 
composited topological excitations in symmetry broken 
states, rotationally invariant spinor BECs and fraction- 
alization of topological excitations in symmetry restored 
states. 

In this section, I will explore the conservation of topo- 
logical charges. Particularly, I want to identify hidden 
topological order in spinor BECs. Furthermore, I gener- 
alize these ideas to the rotationally invariant spin liquids 
and show that topological order coexists with a short- 
range spin correlation, as a result of the topological or- 
der from disorder phenomena. Classification of electron 
liquids using quantum orders was recently reviewed by 
Wen [ |78| . Similar ideas have been implemented for spin 
triplet superconducting liquids WTl. 



A. Hidden topological order in Pontryagin fields in 
BECs 



As discussed at the beginning of the section IV, un- 
der the assumption that the Z2-strings are gapped, one 
can introduce a connection field to characterize all con- 
figurations of n(see also Q,^). The connection field is 
defined as 



A(r) - e^'^e^Vep, 



(50) 



here and n form a local triad; e^-e^ = Sap, ea^e^e^ — 
2n and • n = 0. The field strength is a generalized 
gauge invariant Pontryagin density 



(51) 



dt 



2(Ati - M2) 



(49) 



In (3 -I- l)d, Fq,;, i — x,y,x corresponds to an electric 
field and e^^'^Fij = Hfc represents three components of a 
magnetic field; in 2d, Foi, i = x, y is an electric field with 
two components and the magnetic field has only the z 
component = F^y. 



19 




Fig. 8 The local triad defined at each point r. 

One will be interested in the fluctuations of a topolog- 
ical charge defined as 



Cc;(i)({n(r)}) =^ / dSH.e„, 



47r L 



(52) 



where the integral is carried over the boundary of a large 
volume V and e„ is a unit vector normal to the bound- 
ary. C(7(i) counts the number of hedgehogs in the volume 
V. When the topological charge is a conserved quantity, 
topological different configurations are well-defined dis- 
tinguishable ones. Otherwise, topological distinct config- 
urations are converted into each other by certain instan- 
tons. 



1. Topological order 



I am going to show that at large V, t limit, 
^ hm^ < [Cun) ({n(r)}, t) - Cu^ ({n(r)}, O)] ' >= 



(53) 



in both spin ordered and spin weakly disordered BECs, as 
results of topological long-range order. The topological 
order occurs when hedgehogs are confined, either because 
of the condensation of atoms or an order from disorder 
mechanism. However, in the strongly disordered limit. 



hm < [CuiD ({n(r)}, t) - Cun) ({n(r)}, O)] >= 

t.,V^OC 



tv_ 

Tm 

(54) 



where is finite depending on the tunneling rate of mag- 
netic monopoles(see below). The averages in Eqs,53 5^ 
are taken over the many-spin wave function ^'({n}). In 
2d, we define 



C^(i)({n(r)}) J dxdyU,, 



(55) 



More formally, we can introduce an order parameter 
of the Wilson loop type to characterize the hidden topo- 
logical order. The Wilson-loop integral defined as 



U{1) 



(■Pexp {i j> A- dv)) 



(56) 



has different asymtotical behaviors in the large loop limit 
in the presence or absence of the topological order. When 
the topological charge C[/(i) is conserved, the instantons 
which connect topological different configurations are for- 
bidden; the exponent in the Wilson loop integral is a lin- 
ear function of the perimeter of the loop. On the other 
hand when C'j/(i) is not conserved, the instantons are al- 
lowed and one can confirm the exponent of the Wilson 
loop integral is a linear function of the area enclosed by 
the loop. Though the Wilson loop integral was origi- 
nally introduced for pure gauge fields, here we employ 
it for the characterization of the topological order in the 
spinor BEG where gauge fields are coupled with coherent 
matter fields. 

The topological conservation law has an important 
consequence on collective excitation spectra in spin disor- 
dered BECs. One can show that when the instantons are 
present, the connection field has a massive longitudinal 
component (with a mass rric) 



(57) 



And when the topological charge is conserved, the con- 
nection fields have only transverse components. 



k2 



The topological rigidity can be introduced as 



Pt.o. 



(58) 



(59) 



When topological order is present, pt.o. ^ 0; otherwise, 
Pt.o. vanishes, indicating screening of a small connection 
field (electric component). 





(a) 



fb) 



the total number of Skyrmions living in the 2d space and 
Eqs.^,^ hold in the spin ordered and spin disordered 
BECs respectively (V is replaced by an area S). 



Fig. 9 The responses to an external connection field in 
the absence a) or presence b) of topological order. 
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2. Topological order and the spin stiffness 

To study the conservation of the topological charges, 
we first obtain the equation of motion for using 
Eqs.|,0 



dt 



(60) 



Eq.^ demonstrates that the topological charge is con- 
served if singular space-time events are not allowed so 
that the product n • dn in the right hand side of the 
equation vanishes identically. 

And the topological charges fluctuate only when singu- 
lar space-time events are permitted in the ground state. 
This is rather explicit in 2c?. First, we notice that in 
space X = (i,r 
Gauss equation 



— e'^^Tpi, is the solution of the 



9^H^(x) =^g„,<5(x-x") 



(61) 



in the presence of space-time monopoles {x™} in {2 + l)d 
Euclidean space. 

The space-time monopoles represent the instan- 
tons, which essentially connect a trivial vacuum to a 
Skyrmion configuration. In Euclidean space, a space- 
time monopole located at a; = (to,ro) terminates a 
Skyrmion centered at Tq at time to- There are two im- 
portant aspects in Eq.^. The first is that the topo- 
logical density is conserved in the absence of the space- 
time monopoles in (2-1- 1)D. Second, each monopole event 
causes a change in the topological charge C(7(i) by one 
unit; indeed, following Eqs.|55|,|6l|, 



dr 



(62) 



where the surface contribution has been neglected since 
we will be interested in the leading contribution to Cj/(i) 
per unit square. 

1 1 t t tl t t 1 1 t t 



fb) 



Fig. 10 The spin configurations at a) i = -l-oo, and b) 
t = —oo boundaries of a space-time monopole event; the 
orientation of a spin is indicated by an arrow. 

The probability of finding deconfined space-time 
monopoles depends on the energy of the skyrmions of 



polar coordinates (r, 0p), a static Skyrmion is a configu- 
ration with 

n{r,(j)p) = (sin0(r) cos 0p, sin0(r) sinc/jp, cos(r)) ; (63) 

6{r) varies from at r = to tt at r ^ ^, with ^ an 
arbitrary parameter. For spin ordered BECs where the 
spin Josephson effects should be observed, the energy of 
a Skyrmion is proportional to SttTi^ p/2M and is scale in- 
variant, as reflected in Eq.^ as well. Without lossing 
generality we consider a Skyrmion of a given size ^ Jzo] . 
The connection field is concentrated in a region of the 
size ^, 



-sm6l r — — ^ 
r or 



(64) 



and gets screened at a length scale larger than ^. In 
spin ordered BECs, a Skyrmion is nondegenerate with 
respect to a trivial vacuum. As one will see, this leads to 
the confinement of space-time monopoles. 

To illustrate the point of the confinement, one consid- 
ers ^[/(i) as a function of time measured in some dis- 
crete units and find the following binary representation 
for monopoles. 



a) ...00000000000000000.. 

b) ...00000001111111111.. 

c) ...00000111111100000.. 



U{1) 



1 with respect to a ^(7(1) = configuration. In 



Fig. 11 In a binary representation, monopole- like instan- 
tons are represented by kinks living on a string with 
C[/(i) , the topological charge as the order parameter. 

Bits labeled as 0, 1 as shown in Fig. 11 represent the 
topological charge read out at each moment; a 1-bit cor- 
responds to a Skyrmion while a 0-bit is for a Ct/(i) — 
configuration. The a) string represents a trivial vacuum 
where C'y(i) = at all time. In the second binary string 
b) , a 1-string stands for any scale invariant Skyrmion cre- 
ated at the time to defined by the domain wall separating 
the 0-string and 1-string. So a domain wall in our binary 
representation stands for a monopole Qm = 1, terminat- 
ing a Skyrmion configuration (C(7(i) = 1) at the interface 
and increasing C[/(i) by one unit. A domain wall-anti 
domain wall pair in the third line c) corresponds to a 
monopole- antimonopole pair with Qm = ±1; a 1-string 
here, or a Skyrmion, is terminated at a monopole at one 
end tl and at an antimonopole at the other end 

Since the energy of a Skyrmion is higher than that of 
a trivial configuration, the " energy" of a 1-bit is positive 
with respect to a 0-bit. The energy of a domain wall in 
Fig. 11 which is the number of 1-bits in the structure, de- 
fines the action of the monopole event and is proportional 
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to the length of the 1-string. The action of an isolated 
space-time monopole event, or the creation of a Skyrmion 
therefore is proportional Lthp/2AI, being infinity. {Lt is 
the perimeter along a temporal direction.) The action 
of having monopole-anti-monopole pairs is proportional 
to (<i — 12), the time interval between the creation and 
annihilation of the skyrmions, thus the monopoles are 
confined. For this reason, the right hand side of Eq.|6^ 
vanishes and topological order as shown in Eq.^ is es- 
tablished. 

The conclusion arrived so far can be extended to 3d 
straightforwardly. The change of the topological den- 
sity is caused by the quantum nucleation of monopoles 
instead of Skyrmions. In the binary representation in 
Fig. 11, a 1-bit stands for a static monopole and a 0-bit for 
a trivial configuration. A domain-wall represents termi- 
nation of a static monopole at certain time. Since the en- 
ergy of a monopole is proportional to the system size L in 
a spin ordered BEC, the monopole is nondegenerate with 
respect to the trivial vacuum. Following the same argu- 
ment carried out in 2d, the action to have the monopole 
nucleated is proportional to LLttip/2M , and is infinity. 
This result has an interest impact on the monopole con- 
finement. The interaction between a monopole and an 
antimonopole has to be linear in terms of distance be- 
tween them. The action of a monopole pair nucleated at 
ti and annihilated at t2 with the largest spatial separa- 
tion as |ri— is proportional to |ri— r2|(ti — 12) and the 
monopoles are confined due to the spin stiffness. There- 
fore the right hand side of Eq.^ has to be zero because 
of the confinement, and the topological charge C'[/(i) is a 
conserved quantity. We once again arrive at Eqjs^ in the 
spin ordered BEC. 



3. Order from disorder 

The situation in the disordered limit is more delicate 
and depends on dimensionalities. In the spin disordered 
BEC, the spin stiffness is renormalized to zero and the 
spin fluctuations are gapped. However, the spin fluctu- 
ations do depend on the background connection fields. 
This is explicit in Eq. B3 in the appendix B. Upon in- 
tegrating over spin wave excitations, based on a general 
consideration of the gauge invariance and the parity in- 
variance, we conclude the NLcrM should be reduced to 



1 
2ii 



(65) 



(71 (As) is a function of As, the spin gap measured in units 
of \/ C2jO^/-^/2M^ M]. And in a large-N approximation. 
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2d. 



(66) 



We will analysis the topological order based on Sg. The 
energy of a C[/(i) -configuration is a{L)C^^^^ if the quan- 
tum tunneling is neglected, with a{L) being a function 
of the system size. 

In the absence of the spin stiffness, the energy of a 
Skyrmion or a monopole is determined by its interac- 
tion with the spin fluctuations. This contribution to the 
energy of a skyrmion excitation in a Neel ordered antifer- 
romagnet was previously considered by Auerbach et al. 

. The action of topological instantons was estimated 
by Murthy and Sachdev, Read for quantum disordered 
antiferromagnetic spin liquids [^2|j83| ]. 

In 2d, with the induced interactions, the Skyrmion en- 
ergy is no longer scale invariant and is minimized at the 
size ^ = 00. The connection field of a Skyrmion spreads 
over the whole 2d sheet. The energy of a Skyrmion scales 
as and vanishes as the system size L goes to infinity. 
This implies that the skyrmion configurations become 
degenerate with a trivial vacuum and the monopoles are 
deconfined. This mechanism of deconfinement resembles 
the liberation of fractionalized quasiparticles in one di- 
mension polymers . In one dimensional polyacetylene, 
the ground state has a two- fold degeneracy because of the 
Peierls instability and the domain wall solitons become 
free excitations. 

Unfortunately, when the instantons are liberated and 
the topological charge C[/(i) is not conserved, the en- 
ergy of a configuration Cj/(i) is ill-defined. A more seri- 
ous consideration involves the evaluation of the action of 
monopole-like instantons. In (2-f l)d, the partition func- 
tion of a monopole configuration {x™} represents point- 
like "charges" interacting via long-range Coulomb inter- 
actions. The result suggests that the skyrmions always 
condense in the spin disordered BEC and the topological 
charge Cu(^i^ is not a conserved quantity. It is impor- 
tant to realize that positive and negative Skyrmions can 
be physically and homotopically distinguished because 
of the coherence of the BEC, unlike the situation in a 
classical nematic liquid crystal. 

By moving a Skyrmion around a 7r-disclination, the 
condensate acquires a 7r-phase with respect to the origi- 
nal one. The 7r-phase difference, which can manifest it- 
self in a Josephson type of effect, is one of the signatures 
left behind by the positive or negative Skyrmions. One 
can also distinguish ± Skyrmions by looking at the local 
connection field. In a positive Skyrmion configuration, a 
spin-^ collective excitation, which carries a half charge 
with respect to the connection fields, experiences a con- 
nection field of an opposite sign compared to that of a 
negative Skyrmion This is similar to the situation 
of a point defect discussed in section III. 

So one is able to show that the change of C(7(i), dtC, 
has a short-range temporal correlation because of the in- 
stanton effects. 
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< dtCu(i){t)dtCuii){()) >= exp(- — ). (67) 

And is proportional to L'^^~^t~^ exp{—So)- Imme- 
diately, one recognizes that Eq.^ simply indicates a ran- 
dom walk of ^[/(i) as a function of time in spin disordered 
2d BECs. Therefore , one concludes EqjsJ holds. 

But in 3d weakly disordered BECs where the stiffness 
has vanished, a static monopole with Hfc = r/|r|^ still 
carries a finite energy because of the spin fluctuations in- 
duced interactions. In 3d, as illustrated in 5s(F^i,), the 
spin fluctuations discriminate topologically different con- 
figurations. Particularly, the fluctuations are strongest 
in a Cij(i^ — configuration so that the energy is low- 
est in the topologically trivial configuration. The energy 
of a monopole ^[/(i) = 1 is finite with respect to the 
C'c/(i) = configuration. 

The situation differs from the 2d spin disordered limit 
where the energy of a topologically nontrivial configura- 
tion with C[/(i) = 1 is the same as that of Cc/(i) = one. 
For a monopole, each spherical shell of radius R can be 
viewed as a Skyrmion squeezed into a 2d sheet of size R. 
While the energy of a shell with a large radius is unsub- 
stantial, the smaller shells, which correspond to textures 
of finite sizes and have finite energies, dominate the en- 
ergy of a monopole in the spin disordered limit. Because 
of the singular structure of the monopole, the coupling 
between a monopole and the spin fluctuations lifts the 
degeneracy between the monopole configuration and a 
trivial vacuum in the absence of the spin stiffness. It 
also appears to us that in 3d the spin fluctuations should 
order the topologically different configurations energet- 
ically, independent of the specific form of Ss{F^i,) we 
introduced. The energy of a monopole here is inversely 

1/2 

proportional to gi, i.e. (^EoEt) \nAj^. The action of 
having a monopole, which leads to a change in Cm by one 
unit, is proportional to Lt and is infinite for this reason. 




1 2 3 

Fig. 12 The energy of the monopoles of charge C'f/ji) in 
the absence of spin stiffness. Curve 1 is the result in a 
mean field approximation and curve 2 schematically rep- 
resents the interaction energy between a monopole and 
spin fluctuations in a flnite system. 

The monopoles remain confined even after the S"^- 



rotation symmetry is restored, with the action of having 
a pair of monopoles proportional to {ti — ^2)- Only the 
pair production which conserves Cmi) is allowed. The 
topological order thus coexists with a short-range spin 
correlation in 3d. This can be considered as a case of 
order from disorder phenomena. Though topological or- 
dering occurs in symmetry unbroken states, it mimics 
the conventional order from disorder phenomena as far 
as the role of the fiuctuations is concerned | |85| , |86| , p7t . 
For instance, the spectra of the spin fiuctuations are af- 
fected by a topologically nontrivial back ground in our 
case and by the symmetry breaking in |^,^,^ . Topo- 
logical order is established as a result of the confinement; 
the confinement of the monopoles in this case is driven 
purely by the spin fluctuations, instead of the spin stiff- 
ness in spin ordered BECs. However, by increasing the 
spin disorder further, the degeneracy between C'[/(i) = 1 
and C[/(i) = could be established. This signiflcs the 
deconflncmcnt of the monopoles and the breakdown of 
the topological order. C(7(i) again will have a Brownian 
motion as the time evolves, implied by Eq.|67[ 

B. Topological order in Z2 fields 

A discussion of the hidden topological order of the Z2 
fields can be carried out parallel to the previous discus- 
sion of U{1) fields. For simplicity, I present the result in 
2d case, though a generalization is possible. It is partic- 
ularly convenient to look at the Wilson loop integral of 
the Z2 fields, which takes a form 

Wz. = (n^-'> (68) 
c 

with the product carried out along the loop C. Once 
again, the exponent is a linear function of the perimeter 
when the Z2 vortices are gapped and is proportional to 
the area when the Z2 vortices condense. Consequently, 
the topological charge m 2d defined as 

CzM<^^A)^l[''^, (69) 
c 

is a conserved quantity when the Z2 vortices are forbid- 
den dynamically and is nonconserved when the Z2 vor- 
tices are allowed. 

In the presence of topological order 

lim < [Cz, ({fi,,}, t) - Cz, ({a,,}, 0)] ' >= 0, (70) 

V — *oo 

and in a topologically disordered phase 

9 tv 

lim <[CzAW^l},t)~CzAWv},0)] >=—■ (71) 

In BECs where either the phase or spin is ordered, 
the Z2 vortices are bound with either the half-vortices 
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or the TT-spin disclinations. The energy of the composite 
and the Zi vortices is determined by the stiffness, of the 
phase or spin. The energy of the Zi vortices is therefor 
logarithmic divergent. In the same spirit of the previous 
section, we conclude that the vortices are forbidden 
and confined. The Z^ charge is a conserved quan- 
tity in the spin or phase ordered BECs, as a result of the 
topological confinement. Following this, the Wilson loop 
decays in a fashion of a perimeter law. 

It is perhaps more interesting to look at the possible 
topological order in spin disordered BECs. The spin fluc- 
tuations discriminate the configurations with different Z^ 
charges and the fluctuations are maximized in a triv- 
ial configuration. Therefore, the energy of a Zi charge 
(7^2 = 1 is higher than a trivial one of Cz^ = 0, though 
the spin stiffness vanishes. 

Indeed, because of a coupling between the Z2 fields 
and the spin fluctuations, an effective interaction of the 
Z2 fields is induced upon an integration of spin fluctua- 
tions, 

^---^Xli^rr'; (72) 

f/2 is a function of the spin gap developed in disordered 
limit Eqj7| once more illustrates a general idea of 

the topological order from disorder. A Z2 vortex Y[o '^ij 
does acquire an energy, by interacting with the spin fluc- 
tuations. The energy of a Z2 vortex in spin disordered 
BECs depends on the spin gap; it is finite but enormous 
in the spin weakly disordered limit in 2d, reflecting an 
energetic difficulty to have Z2 vortices due to the weakly 
gapped spin fluctuations. The rest of discussions can be 
carried out in a binary representation. The action of 
an 1-string, representing termination of a Z2 vortex is 
infinity. The instantons connecting configurations with 
different Z2 charges are suppressed and is conserved. 
Topological order can be established by the spin fluctu- 
ations. 



VII. SYMMETRY RESTORING IN FINITE SIZE 
SPINOR BECS 

Since the experiment was done in BECs of a few 
millions "^^Na atoms, it is also particularly interesting 
to consider the symmetry restoring due to a finite size 
effect. One takes a weakly interacting limit where quan- 
tum fluctuations of finite wave lengths are negligible. In 
a zero mode approximation, one can neglect the spatial 
fluctuations of n and obtain from Eq.|3^ 

^ , , ^ , , 4:tc2ph^ f . sin^ X 
C^(q = 0, i) - C^(q = 0, 0) = ^ J dx—^ 

(73) 



which is valid when the result is less than unity. C,,(q = 
0, 0) should be considered of order unity. Formally speak- 
ing, the extrapolation of the above result to t — > 00 shows 
that the correlator diverges as t goes to infinity. This im- 
plies that n{t), at i ~ N/c2p{set by the two-body scat- 
tering length), becomes uncorrelated with the original 
orientation of n(0). At longer time scales, n precesses 
an extremely slow quantum diffusive motion on the unit 
sphere. That is to say, the nematic order parameter can 
rotate freely in the parameter space, as seen from C,,(0, 0) 
as well and restore the broken rotation symmetry when 
N is finite. 

Following Eq.|l^, in the OD BEC the energy gap of the 
low lying excitations scales of at large iV limit. In 
the ground state, < n because of the rotational in- 
variance. Consider a wave packet with n confined within 
a region centered at no = e^. The spread V < 6'^n >o 
on the unit sphere at t = is very small. The spread 

< (5^n > is defined as the expectation value of + iiy- 
In spherical polar coordinates {9, (p) defined on a unit 
sphere, L is a differential operator in terms of {9, (p) and 
the eigenstates of are spherical harmonics Yi^m{9, 0) . 
The corresponding wave packet can be constructed as 

l=2n 

B = ^exp(-i^),<52„>o=^ (74) 
I 

where constant Aq is estimated in the Appendix C. 
^i,o(^i0) are spheric harmonics. Only states with even 
number / are present in the wave packet (we assume N 
is an even number at the moment). For this packet, 

< rij; >=< rij, >= and < >w 1 when ct ^ 1 
but a /N is vanishingly small. 

The energy and the spread in L can be estimated as 

AE='^-^^,<5^l.> = AA,a (75) 

according to the derivation in the Appendix C. There- 
fore, an excited state where n is well defined on the unit 
sphere up to a spread V < >o has an infinitesimal 
energy cost when N approaches infinity. 

The wave packet constructed in Eq.|7j, being energeti- 
cally extremely close to the singlet ground state, indicates 
that the singlet is very unstable with respect to external 
magnetic fields. A field of the order o(l/A^) will stablize 
the symmetry broken state in Eq.|^. However, when one 
deals with BEC of a small number of atoms, the level 
spacing in the excitation spectrum gets bigger. The sys- 
tem becomes perturbative with respective to an external 
magnetic field, as far as the disturbance is less than the 
energy gap in the problem. So the claimed ultrasensitiv- 
ity is simply a consequence of nearly degenerate low lying 
excitations in the large N limit. Moreover, this ultrasen- 
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sitivity is absent in the spin disordered states where a 
many body spin gap is opened up. 

According to the equation of motion in Eq.|ll|, as the 
time evolves, because of the two-body scatterings with 
C2 > 0, n{t) gradually precesses along the finite SL. In- 
deed, the spread of n at a time t becomes 



(76) 



which is valid at t < Tc = a/AoAE. In an opposite limit 
the spread is of order unity. 

Following EqJ?^, at a time of order 1/AE, \J < S'^n >t 
exceeds the initial spread < S'^n >o- At a longer time 
Tc, n starts to explore the whole unit sphere S'^ and the 
rotation symmetry is restored quantum mechanically in 
the presence of spin-dependent scatterings (see Fig. 13). 
A discussion is present in appendix D when an external 
field is applied. 





(a) 



(b) 



(c) 



Fig. 13 The symmetry restoring caused by the spin- 
dependent interactions (c2 > 0). Schematically shown arc 
the spreads of n on a unit sphere at t — (a), 1/AE (b), 
1/AE X 1/V< S'^n >o(c). 



Eq.^ also imposes restrictions on the measurement of 
a symmetry broken state. An measurement of n in BECs 
of "^^Na discussed here excites the ground state to an ex- 
cited state where n has a finite spread on S^; the energy 
of this state is infinitesimal small in the thermal dynam- 
ical limit. The time scale for n to relax, or to restore the 
broken symmetry is determined by the two-body spin de- 
pendent scattering lengths and the number of atoms in 
BECs following Eq.^. Experimentally, the relevance of 
the singlet ground state or the symmetry broken state 
depends on the life time of the metastable gas or the in- 
terval between measurements. Any measurement taken 
at a time scale longer than (t/AqAE should reveal the 
symmetry restoring because of zero point motions of n. 
Possible measurements, including measurements of bire- 
fringence and small angle light scatterings, to detect or- 
ders in n in the weakly interacting limit were discussed 
in |7|. 

It is worth emphasising that by contrast to a classi- 
cal nematic liquid crystal where n is pinned along the 
external field and all the Goldstone modes are gapped, 
n can still freely rotate in a plane perpendicular to a 
weak magnetic field in the quantum spin nematic state 
under consideration. So the singularity in small angle 
light scattering survives an external magnetic field; only 
the amplitude is reduced by a factor of one-half, when an 
external field is applied. In a nematic disordered phase, 
the singularity gets smeared at q ^ is the correla- 

tion length) , which can be used to study a possible phase 
transition. 



Eq.^ shows how a symmetry broken state in finite 
size BECs of '^^Na becomes a rotationally invariant state 
at an asymptotically long time. For ^^Na atoms inter- 
acting via a two-body interaction with C2 > 0, it was 
argued in Ref. [^2| that in a zero mode approximation 
the atoms should condensate at \l,Tn = >, with the 
rotation symmetry broken. On the other hand, based on 
an analogy to the nonlinear wave-mixing theory, it was 
shown in Ref. |Q that the exact ground state should be a 
spin singlet with the rotation symmetry preserved. The 
symmetry restoring in a finite size system as discussed 
above explains the apparent inconsistency between these 
previous theoretical works. 

Given the number of particles as 6 x 10^ and /Iq ~ 
100Hz{5Q0nK), the typical time scale for this to take 
place is tji — NIq/Tt? ~ Iday. This time scale is much 
longer than the life time of the trap Q. Therefore the 
features associated with the singlet wave function pro- 
posed in Ref. |3^ can not be observed under that partic- 
ular experiment condition. However, for a smaller cloud 
with 10'^~* atoms, the symmetry restoring can take place 
within 10— lOOsecs, before the recombination takes place. 
In this case, the ground state should be characterized by 
a singlet wave function, instead of a saddle point solution. 



VIII. DISCUSSION 



A. Relation to an HAFS 



Though the NLcrM description of the quantum spin 
nematic states bears some similarities with that of the 
Heisenberg antiferromagnetic spin system(HAFS), there 
are at least two differences between the present problem 
and the Heisenberg problem. First, for the NLcrM of the 
HAFS, the ground state is invariant under a rotation on a 
sphere the order parameter space is a unit sphere . 
iTi{S'^) is homotopic to zero and the linear singularity is 
absent in this system. 7r2(S'^) — Z and the corresponding 
point defects are usual hedgehogs with a winding num- 
ber N distinguishable from ~N . Overall, the topological 
defects in the Heisenberg spin system are qualitatively 
different from what we have in the QSNS described in 
section HI. 

However, when a weak magnetic field is applied, n in 
the NLcrM of the HFAS lives on a unit circle and 
7ri(S'^) = Z i& an integer group. The linear defects are 
integer vortices. Consequently, in 2d in the presence of an 
external field, the interaction between topological linear 
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defects in the QSNS at an asymptotically large distance 
decays in a same fashion as that in the HAFS. We expect 
that in this specific limit, the finite temperature phase 
transitions are of the same universality classes. In addi- 
tion, TT-spin disclinations superimposed with half vortices 
interact with each other via one half of the interaction 
strength between two integer vortices. Because of the 
reduced interaction strength between half vortices in the 
QSNS at H 0, the critical temperature is lower than 
that of the HAFS. 

It should be noticed that though the spin-dependent 
scatterings lift the degeneracy between different spin 
configurations and result in a spin nematic state when 
C2 > 0, the critical temperature is determined by the 
density, i.e. Tc ~ h^p'^^^/M which is independent of C2 
in the leading order approximation. When the zero point 
motions of n in weakly interacting BECs are taken into 
account, Tc is shifted downward by an amount o(/) pro- 
portional to C2. 

Second, in the Heisenberg antiferromagnetic system 
with a underlying spin S, both the zero point kinetic en- 
ergy (JS) and the potential energy ( JS*^) depend on the 
exchange integrals and the corresponding /, the dimen- 
sionless coupling constant in the NLcrM, is inversely pro- 
portional to the spin S only. This makes observing zero 
temperature order-disorder phase transitions extremely 
difficult. In the current situation, the potential energy 
is independent of the interaction while the zero point ki- 
netic energy is proportional to the scattering length. / 
is determined by the density and the scattering length 
instead of spins of Bosons and is a continuous variable in 
the QSNS. 

B. Classical nematic liquids 

The classical nematic order-disordered transitions were 
studied a few years ago based on a lattice gauge theory 
[^3|. The inversion symmetry of a nematic director is 
incorporated as a Z2 gauge invariance in a lattice model, 

— = -J^arr'nr ■ Ur' - K^Y\_'^rr' (77) 
II ^j./ □ □ 

where arr' — il are introduced as the Z2 gauge fields. 
The sum is again over all elementary plaquettes. The 
Hamiltonian is invariant under a Z2 gauge transforma- 
tion: rij. — !■ — rij,, a^r' — —Urr' (all rr' links involving 

r). _ 

The Hamiltonian in Eq. ^ depends on two energies 
scales: J as the microscopic nematic interaction and K 
as the core energy of a tt spin disclination. At K ^ 
the Z2 field is activated. As J decreases, a spin nematic 
state is driven into an isotropic state where the free en- 
ergy of an infinity long tt disclination is finite. The phase 
transition is of the first order and believed to be the one 



described by the Landau theory. Ki K ^ 00, the Z2 
gauge field is frozen. In this case, as J decreases, the ne- 
matic state is driven into another isotropic state where 
the free energy of a tt disclination is proportional to the 
length of it. The transition is of the second order. 

The gauge theory description of the classical nematic 
liquid predicts two isotropic and optically transparent liq- 
uids. They are different only in topological order. One of 
them has the hidden Z2 order and the other one doesn't. 
A recent interesting experiment by Tokyo liquid crystal 
group did support this scenario |Q. The optical tur- 
bidity measurement showed a transition between an op- 
tically turbid phase, which apparently can be identified 
as a nematic phase and an optically transparent state 
which represents an isotropic liquid. However, within the 
isotropic liquid, the dynamical slow down measurement 
indicates an additional phase transition between two op- 
tically transparent phases. Though it is not clear at the 
moment whether these two phases carry different topo- 
logical order, the situation there is quite encouraging. 

C. Fermionic nematic liquids 

Finally, I want to remark on other electronic liquids 
which are believed to be "nematic" in nature. The first 
one is the ground state of Sr2Ru04, which was recently 
confirmed to be a p-wave superconducting state in a se- 
ries of beautiful experiments |]89| , po| , ^ , p2| . In a p-wave 
superconductor of Sr2RuOA, a Cooper pair is a spin 
triplet and experiments further suggest that it be an 
Anderson-Morel state If a pair is considered as a 

composited boson, then it occupies a to^^^ = channel of 
spin one states, up to a rotation of the quantization axis. 
The order parameter can be written in a tensor form as 

Aap = Ao(T) exp(ix) (^2(1 • n)o,p{\<ix ± i^y) (78) 

in the absence of spin-orbit couplings, with Ui, i — 0,1,2 
being Pauli matrices. By examining this pairing wave 
function, one expects that many aspects of Sr2RuOi 
should be very similar to the spinor BEC described 
above. In particular, the concepts of the hidden sym- 
metries and hidden gauge fields should be applicable in 
those systems Of course, in addition to this, there is 
a new twist in p-wave superconductors, that is the broken 
time invariance and broken parity p^ . It leads to pos- 
sible statistical transmutation of topological excitations 
and quasiparticles. This remains to be understood in the 
context of quantum nematic states discussed in this arti- 
cle. Results on spinor BECs discussed in this article will 
have no doubt to be demonstrative for the understanding 
of the physics in Sr2RuOi superconducting crystals. 

The Fermionic nematic liquids were also suggested for 
cuprates recently. The nematic model discussed in [ [73[ , 
and the QSNS in the spinor BEC discussed here are 
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intimately connected with the recent Zi gauge theory 
developed for high- Tc superconductors |Q. To under- 
stand the connection between a " fractionalized" metal 
(non Fermi liquid) and a normal metal(Fermi liquid), the 
authors of that work pointed out that fractionalization 
and confining issues in cuprates can be partially under- 
stood by taking into account an inversion symmetry in 
a spinon-chargon model. In anisotropic phases of some 
strongly correlated electron systems such as "quantum 
Hall nematic phases", fermionic nematic liquids, or ne- 
matic metals were also investigated by Kivelson, Fradkin 
and their collaborators |94 9^. V. Organnesyan et al. re- 
cently showed the breakdown of the fermi liquid theory 
in electronic nematic phases. 

The investigations carried out for the spinor BEC 
therefore have possible impacts on the understanding of 
spin-charge separation or electron fractionalization phe- 
nomenon, which was first envisioned for cuprates about 
15 years ago [^,^,^,^. As one of the most interesting 
inventions by theorists in condensed matter physics, spin- 
charge separation still continuously receives a lot of at- 
tentions from condensed matter physicists. Though mi- 
croscopically BECs studied here are very different from 
those strongly correlated electronic systems, the long 
wave length physics do bear many similarities to some 
models employed for electron liquids. Studies of the 
quantum spin nematic states, in particular the quantum 
orders and spin fractionalization , not only are im- 

portant for the understanding of spinor BECs, but also 
shed some light on some long standing issues in many- 
body physics. Being free of imperfections and easy to be 
manipulated magnetically or optically, the QSNS might 
be a potentially new platform for exploring some basic 
ideas in strongly correlated electron systems, which have 
not been settled. 



APPENDIX A: SPIN WAVE FUNCTIONS IN 
{V,V) REPRESENTATION 

In (m, v) representation, states of spin S are polynomi- 
als of degree 25* in terms of u, v. Examples for S* = 0, 1, 2 
are given below 

15 = 0,52 = >= w" ^ 1; 

\S^-,S. = - >= ^-^.^ ^-cos-exp(z-), 

= 2'^^ = —2 >= V r ^ V 2 2 '^P(-*2)' 



IS" = 1, S'^ = 1 >= \/3ii^ ^ VScos^ - exp(i(/.), 
\S — \^Sz^^>= VGuv a/6 sin 0, 

\S = 1,S, = -1 >= -VSv'^ -VScos^ - exp(-i(/)). (Al) 



Finally, 



V6u(n)w(n) = cos 6*11,0 > + 

— = sin0exp(— z0)|l, 1 > + 
v2 

— r= sin0exp(z(/))|l, — 1 > 
v2 



(A2) 



where n = {0,(p). At 9 — tt/2, the amplitude of 
\S — 1, Sz — > state vanishes. 



APPENDIX B: DERIVATION OF THE RG 
EQUATION 



Following [p8[ , we introduce the unit vector n in a form 



of 



n = (1 - ^ 0')'/'eo + J2 '^'J^r,, 
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(Bl) 



7y = 1, 2. <f>n can be considered as spin waves (fast vari- 
ables) around given slowly varying 69(2:). 

Given the constraints on eo,e^ in Eq. Bl, we further 
introduce 

a^eo = SO-'e,,, d^B^ = A^B^, - B^^eo. (B2) 

Inserting Eqs.Bl, B2 into Eq.^ and expanding (1 — 
(jy^) in terms of 0^,, one obtains 

1717' n 

(B3) 
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By integrating out fluctuations of 4> from e~' < |p| < 1 
and e~' < |e| < 1 and restoring the momentum and the 
energy range, we obtain recursion relations for f(J,) in 
d + 1 dimension 



f{l) 

n(0 



n(0], 



7(0) 

^ d'^q da; 1 



(B4) 



(B5) 



{271^ 27r ^2 + q2- 
The RG equation foUowing Eq.B4 is 

^ = (l-d)/ + i?,/2, 

with Ri = l/27r,i?2 = l/27r2,i?3 = 1/4tt'^ 



APPENDIX C: CALCULATIONS OF SPREADS 
AT A ZERO FIELD LIMIT 



In spherical coordinates, 

d d 
La; = ih{sm (j)— + cot e cos (/>7r-), 
o& dm 



d 



d 



Jjy — ih{— COS 0— + cot 6 sin 4)^;-) 
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(CI) 



Taking into account Eq.Cl, we obtain 
<{5'n) >,= 1^(1-Mu0x 



B 

I 

exp it exp 

4ct 2iV ' 



/2 



4ct 



it 



(^ + 1)' 



2A^ 

^2 



(2/ + l)(2Z + 3) (2/ + l)(2/-l)' 



. (/ + !)(/ + 2) Z(/-l) 



^(2^ + l)(2/ + 3)^ 



-(2; + l)(2/- 1)' 



(C2) 

In the limit ct ^ 1, we obtain results in Eq.^ with 

(C3) 



APPENDIX D: CALCULATIONS OF < (5^N >t AT 

7? / 

At iJ 7^ 0, the equation of motion becomes particu- 
larly simple. L — > = ifid/d9 and n — (cos0, sin0); 6* is 
the coordinate on a unit circle. 



A wave packet of interest where <n2:>~l, <nj,>= 
can be written as 

vl'(n,t) = i exp(-i^)exp(zZ0-ziZ2£2^) 



i=0.±2... 



(Dl) 



where cr 3> 1 is the spread of Ig in the wave packet. One 
can again show that the energy, spin and the correspond- 
ing spread of n are 



N 



< S'^Uy >(= ylo[- +4cr( 



tC2p, 

N ' 



(D2) 



Ao(cr) is a quantity of unity. The last equation in Eq.D2 
is valid when t < l/AE^/a. In deriving Eq.D2, we em- 
ploy Eq.C2 and notice that 



(D3) 



in this case. 
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